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I.  INTRODUCTION 

TRIXY  is  a  computer  program  that  provides  a  multigroup  analysis  of 
nuclear  reactors  in  three  space  dimensions .  The  program  can  be  run  on  an 
IBM  704  computer  that  has  at  least  8192  words  of  memory  and  six  tape  units. 
The  input  data  consist  of  cross-sections,  geometrical  configurations  describing 
the  location  of  material  within  the  reactor,  control  rod  positions,  mesh  spacing 
for  the  three  axes,  and  boundary  conditions.  For  each  energy  group,  g,  the 
program  computes  ^  ,  the  neutron  flux,  by  a  numerical  approximation  to  the 
age-diffusion  equation.  When  has  been  calculated  for  all  groups,  the  fission 
source,  lb  ,  is  computed.  The  criticality  and  source  error  are  calculated 
and  printed.  If  the  fission  source  has  not  converged  to  the  degree  of  accuracy 
requested,  1^  is  recomputed.  The  results  of  the  program  are  edited  so  that 
values  of  flux  and  power  density  for  any  section  of  the  reactor  can  be  printed. 

GENERAL  FEATURES  OF  TRIXY 

1.  Three- dimensioned  Cartesian  coordinates  are  used. 

2.  Variable  mesh  spacings  may  be  chosen  on  all  three  axes. 

3 .  There  may  be  up  to  72  mesh  planes  perpendicular  to  any  of  the  three 


axes . 
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4.  Ur  ’  six  energy  groups  may  be  used, 

5.  Up  to  975  different  material  compositions  are  permitted  within  the 
mesh  (with  32,  768  word  memory). 

6.  The  mesh  is  bounded  by  a  rectangular  parallelepiped.  Boundary  condi¬ 
tions  are  specified  independently  for  each  of  the  six  faces,  either  zero 
flux,  or  zero  or  non- zero  current. 

7 .  The  product  of  the  number  of  mesh  points  and  the  number  of  energy 
groups  is  limited  to  150,  OCX). 

8 .  The  contribution  to  the  neutron  source  due  to  slowing  down  of  neutrons 
is  restricted  to  degradation  from  the  next  hi^er  energy  group. 

9 .  Fluxes  in  each  energy  group  are  calculated  by  an  iterative  technique 
(accelerated  Liebmann). 

10 .  The  edit  routine  yields  fluxes  and  power  densities  for  mesh  points  in 

(3) 

selected  columns  (a  column  is  a  mesh  line  parallel  to  the  *  axis), 
the  power  integral  and  average  power  density  for  selected  columns,  flux 
Integrals  and  average  flux  values  for  each  composition  region  of  the 
mesh,  and  the  number  of  fissions  per  source  neutron  (proportional  to 
criticality) . 

11.  A  power  guess  for  the  first  run  (a  run  is  a  power  iteration)  may  be 
specified  in  a  coarse  mesh,  or  the  results  of  an  earlier  calculation  may 


be  used  if  the  number  of  mesh  planes  in  each  direction  is  identical  for 
the  two  problems . 

K  G 

12.  Adjoint  fluxes  may  be  computed,  provided  the  ratio  F  /F  is 

independent  of  position  throughout  that  portion  of  the  mesh  that  contains 

Q 

fissionable  material.  F  is  the  fission  cross  section  for  energy  group 
g,  and  G  is  the  lowest  energy  (thermal)  group . 

LIMITATIONS  DETERMINED  BY  SIZE  OF  MEMORY 

Although  TRIXY  was  written  primarily  for  use  on  a  32,  768  word  com  ¬ 
puter,  it  may  also  be  run  on  a  8192  or  16, 384  word  computer  with  the  following 
input  restrictions ; 

Words  of  memory  8192  16, 384  32, 768 

Product  of  the  number  of  mesh 

(2) 

planes  perpendicular  to  the  X  and 

X^^^  axes.  1000  2000  4000 

Approximate  number  of  mat¬ 
erial  compositions  permitted  within 

the  file  of  cross  sections  records.  500  1500  3000 

Number  of  material  comp¬ 


ositions  permitted  within  the  mesh. 


200 


475 


975 
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If  there  is  a  change  in  mesh  size  or  material  composition  along  any 
column  at  only  every  fifth  point,  tlien  the  running  times  for  a  given  problem 
will  be  approximately  the  same  for  any  size  machine. 


5 


11.  DESCRIPTION  OF  THE  MESH 

GEOMETRIC  PROPERTIES 


Consider  the  Cartesian  coordinate  axes  and  Mesh 

planes  are  perpendicular  to  these  axes  and  have  coordinates  X^^l  X^^^,  X^j^^ 

1  j  X 

The  integers  i,  j,  and  k  identify  the  planes  and  have  the  maximum  values 
I.  J,  K  =  72. 


Mesh  points,  identified  as  (i,  j,k),  are  located  at  the  intersections  of 
the  mesh  planes  of  corresponding  indexes.  Point  (1, 1, 1)  is  located  at  the 
origin  of  the  coordinate  system. 

The  spacing  (cm)  between  adjacent  mesh  planes  need  not  be  uniform 
and  is  defined  by; 


1 

- 

1+1 

-  x<« 

1 

for 

1  <  i  < 

I  -  1; 

O 

U 

<3 

(Ill) 

a<2) 

J 

-  y<2) 

■  Vi 

-x<^> 

J 

for  1 

<j  <J-i; 

0 

f  = 

0, 

(II  2) 

^k 

(3) 

=  ^k+l 

1 

X 

for  1 

V 

M 

V 

'^0  - 

[> 

X  w 

II 

0. 

(II  3) 

The  six  mesh  planes,  i,  i+1,  j,  j+l,  k,  and  k+1  bound  a  rectangular 
parallelepiped  of  uniform  composition.  This  parallelepiped  is  identified  as 
"composition  box"  (i,  j,k)  and  its  volume  is: 


6 


(II  4) 


A  composition  box  is  divided  into  ei^t  boxes  of  equal  volume,  identified 
as  octants.  The  octant  nearest  mesh  point  (i,  j,k)  is  identified  as  octant 
(i,  j,  k),  and  its  volume  is: 


V. 


i.  j.k 


=  A 


(1)  ^(2) 

i  j 


Af/8. 


(II  5) 


Each  mesh  point  that  does  not  lie  on  a  boundary  plane  is  surrounded  by 
eigjit  octants  which  may  be  of  different  composition.  Boundary  points  touch  one, 
two,  or  four  octants,  respectively,  for  points  that  lie  on  corners,  edges,  or 
faces  of  the  mesh .  The  total  region  of  all  octants  that  touch  a  given  mesh  point 
(i,  j,  k)  is  defined  as  "flux  box"  (i,  j,  k),  and  this  volume  is : 


2 


£ ,  m, n=0 


i-i*.  j-m, k-n’ 


(II  6) 


Some  of  the  octants  in  flux  box  (i,  j,  k)  may  contain  fissionable  material. 


The  beta  function  is  defined  as : 


/3 .  .  ,  =  1  if  V.  .  ,  contains  fuel, 
i,],k  i,j,k 

(U  7) 


/3.  .  .  =  0  if  V.  .  , 

i.J.k  i,],k 


does  not  contain  fuel . 
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Hence,  the  fissionable  volume  in  a  flux  box  (i,  j,  k)  is; 

V  =  >  V  B 

i.  j,k  n^m,'n=0  i-J?  ,  j-m,  k-n' 

and  the  total  core  volume  is  • 

I.  J,  K  .. 

P  <^1,  ,  1. J>k 

i..j,k=l  •” 


(11  8) 


(II  9) 


CROSS  SECTIONS 

Each  material  composition  in  the  mesh  is  identified  by  an  integer  I 

where  Z^=l,  ^3=3,  etc.  To  distinguish  between  compositions  that 

contain  fissionable  material  and  those  that  do  not,  it  is  required  that  1  <  i  < 

-  n  - 

999  for  compositions  that  contain  fuel  and  1001  <  I  <  1999  for  non -fissionable 

-  n  - 

compositions .  G  sets  of  values  D®,  A^,  M®,  and  are  associated  with  each 
in  the  fissionable  range,  where  the  superscript  g  denotes  the  energy  group 

g  a 

I  5  S  D  is  the  diffusion  coefficient  (cm),  is  the  slowing  down 

- 1  g  "Iff 

cross  section  (cm  ),  and  F  is  the  fission  cross  section  (cm  ).  A°  is  the 
removal  cross  section  (cm  which  is  the  sum  of  the  absorption  and  slowing 
down  cross  sections .  G  sets  of  values  D^,  A®,  and  are  associated  with 
each  in  the  non -fissionable  range. 

A  library  of  diffusion  coefficients  and  cross  sections  are  maintained 
on  tape.  Each  file  in  this  library  contains  four  separate  records: 
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(1)  The  identification  number  of  the  file,  the  number  of  energy  groups  G, 

JC®  for  each  group  (X  is  the  fraction  of  the  fission  spectrum  that 
occurs  in  group  g),  and  the  number  of  words  in  each  record  in  the  file, 

(2)  Diffusion  coefficients  and  cross  sections  for  compositions  containing 

fuel  in  the  order  of  increasing  /  , 

n 

(3)  Diffusion  coefficients  and  cross  sections  for  non -fissionable  composit¬ 
ions  in  the  order  of  increasing  and 

(4)  Ratios  of  outgoing  current  to  flux,  Z®,  for  various  reflector  composit¬ 
ions,  each  identified  by  an  index  ^  ,  where  2001  <  i  <  2999. 

n  “  n  ~ 

O' 

G  values  of  Z®  are  listed  for  each  and  compositions  are  listed 
in  the  order  of  increasing  The  Z®  values  are  used  to  specify  non-zero 
current  boundary  conditions . 

The  diffusion  coefficients  and  cross  sections  of  composition  box  (i,  j,k) 

are  identified  by  a  composition  index  Z.  ,  where  Z.  is  one  of  the 

1»  i>  J> 

integers  in  records  two  and  three .  Current  boundary  conditions  for  the  bound¬ 
ary  planes  are  identified  by  the  composition  indexes 

^0  h  ^0 

/  0\  /QX  /  \ 

Uy  ,  and  i/„  ,  respectively  for  the  six  planes,  where  ,  is  identified 

with  an  Z^  in  record  four .  Composition  indexes  Z^  ^  ^  for  each  composition 

box  of  the  mesh,  combined  with  a  file  of  diffusion  coefficients  and  cross  sections, 

provide  a  complete  set  of  cross  section  data  for  the  mesh.  Values  D®  ,, 

i,J,k 
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A®  ,  ,  ,  and  for  all  groups  associated  with  index  I 

i.j,k  i.j.k  i,j,k  i.j.k 

describe  the  nuclear  properties  of  composition  box  (i,  j,  k). 

SPECIFICATION  OF  CONFIGURATIONS  OF  MATERIALS 

A  configuration  is  an  area  of  uniform  material  composition  in  a  plane 
(3) 

normal  to  the  X'  '  axis  (a  k-plane) .  List  A  and  list  B  are  used  to  project 
these  plane  figures  into  three-dimensional  volumes .  In  constructing  each 
k-plane,  list  A  is  examined,  and  the  appropriate  configurations  are  placed  on 
the  plane;  then  list  B  (control-rod  list)  is  examined,  and  the  control  rod  con¬ 
figurations  are  written  over  the  list  A  configurations . 

Each  configuration  is  a  record  in  a  file  of  configurations  preceding  the 
files  of  Cross  Section  Data.  The  maximum  number  of  configurations  is  limited 


only  by  the  physical  length  of  the  tape. 
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m.  INPUT  DATA 


CLASSIFICATION  OF  INPUT 


There  are  two  main  divisions  of  input  for  a  problem.  The  first  set  of 
input  consists  of  a  tape  containing  a  file  of  configurations  and  several  files  of 
Cross  Section  data.  It  is  expected  that  one  tape  can  be  used  for  a  number  of 
problems . 

The  second  set  of  input  consists  of  the  following: 

A.  Input  control  data . 

1.  Miscellaneous  data,  such  as  I,  J,  K,  G,  for  each  group, 


<T 

,  Vp,  gp  (F®  =  0  everywhere  for  g  <  g’),  si 

F  F  F  F  F  F 

mesh  planes  (i  ,  ,  i  ,  j  .  ,  j  ,  k  .  ,  and  k  ) 
nun  max  ^mm  max  mm  max 


SIX 


that  boimd  all  fissionable  material  in  the  mesh,  and 


(P) 


for  each  boundary  plane  of  specified  current -to -flux  ratio. 

The  machine  calculates  I,  J,  K,  and  V„.  These  calculations 

r* 

are  listed  as  input  for  the  purpose  of  checking  the  consis¬ 
tency  of  the  input  data , 

2.  List  A 


3. 


List  B 
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4.  Information  for  calculation  of  mesh  spacings.  Input  data 
list  coordinates  along  each  axis  and  the  number  of  equally 
spaced  intervals  between  adjacent  coordinates 

5 .  Input  control  data  for  the  EDIT  routine .  Output  data  that 
may  be  requested  are; 

a.  The  power  density  for  each  mesh  point,  either  for  the 
entire  mesh  or  for  selected  columns, 

b.  The  average  power  density  and  the  power  integral  for 
a  column,  either  for  all  the  columns  of  the  mesh  or 
for  selected  columns, 

c.  The  flux  at  each  mesh  point  for  selected  energy  groups, 
either  for  the  entire  mesh  or  for  selected  columns, 

d.  The  flux  integral  and  average  flux  for  each  of  the 
different  composition  regions  of  the  mesh  for  each 
group,  and 

e.  The  number  of  fissions  per  source  neutron  (f),  and 
the  contribution  of  each  energy  group  to  f. 
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B.  A  power  guess. 

If  a  power  guess  and  flux  guess  are  available  from  a  pre¬ 
vious  problem,these  may  be  used.  If  no  guess  is  available,  a 
rough  power  guess  must  be  supplied  in  a  coarse  mesh. 


MESH  SPACINGS  AND  COORDINATES 

A  series  of  numbers,  n^^^,  ....  0, 

aabb  I  aab 

xj^^^  ...,  X^^\  0,  n^^^  X^^^  n^^^  X^^\  ....  X^^  ,  0,  is  used  to  specify 

mesh  increments .  Along  the  X^^^  axis,  the  interval  between  coordinate  X^^^ 

a 

and  the  origin  is  divided  by  mesh  planes  into  n^^^ equal  intervals.  The  interval 

a 

between  X^^^^and  X^^^  is  divided  by  mesh  planes  into  equal  intervals,  until 

the  first  zero  is  reached.  The  zero  indicates  the  end  of  the  specifications  of 

mesh  spacings  along  the  X^^^  axis .  A  similar  interpretation  applies  to  the 

(2)  (3) 

specification  of  the  mesh  spacings  along  X'  '  and  x'  '  axes . 
CONFIGURATIONS  WITH  LIST  A  AND  LIST  B 


A  hypothetical  reactor  is  constructed  using  the  configurations  in  Figure  1 
to  demonstrate  the  method  of  specifying  configurations  and  the  projection  of  these 
plane  figures  into  three-dimensional  volume. 

Each  of  the  configurations  is  enclosed  by  a  rectangle.  To  describe  con¬ 
figuration  1,  the  number  of  intervals  along  the  j  axis  is  seven.  Reading  up 
the  first  column,  there  is  one  block  within  the  configuration.  There  are  six 
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Configuration  1 


Configuration  3 


Two  dimensional  representations  of  hypothetical  configurations. 


FIGURE  1 
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blocks  outside  (-).  Therefore,  the  first  column  is  described  by  1,  -6,  100, 
where  100  denotes  the  end  of  a  column  „  This  process  is  repeated  for  the 
remaining  columns,  and  1000  is  written  at  the  end  of  the  configuration. 

Configuration  1  (n  =  7)  is  completely  described  by  the  following 
series  of  numbers : 

1,7,  1,  -6,  100,  2,  -5,  100,  4,  ~3,  100,  5,  -2,  100,  5,  -2,  100,  6,  -1,  100,  7,  1000. 
Similarly  for  configuration  2: 

2,  3,  -1,  1,  -1,  100,  3,  100,  -1,  1.  -1,  1000, 
and  configuration  3: 

3,  1,  1,  1000. 

In  this  hypothetical  reactor  (see  Figure  2)  let  configuration  1  be  the 

core,  composition  index  I  =1,  and  extend  from  plane  k=3  to  k=5  where  the 

n 

configuration  changes  to  composition  index  2.  This  composition  index  extends 
to  k=7.  The  core  is  surrounded  by  composition  =  1001  on  all  sides  except 
i=I  and  j=l. 

A  control  rod,  configuration  2,  composition  index  1002,  extends  from 
k=4  to  k=K.  The  volume  under  this  control  rod  is  composition  =  1001. 

Another  control  rod,  configuration  3,  composition  index  1003,  extends 


from  k=3  to  k=K. 
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List  A  k  C  ?  i  i 

n  n 

1  0*  1001  1  1 

3  113  1 

5  12  3  1 

7  0  1001  1  1 


List  B  k  C  Z  i  i 

n  n 

3  3  1003  9  1 

3  2  1001  6  2 

4  2  1002  6  2 


ROUGH  POWER  GUESS 


Coordinates  (cm)  of  a  coarse  mesh  are  listed  as  a  series  of  numbers 

....  0,  ....  X^p.  0,  xp.  ....  xP ,  0.  The  values 

entered  for  xj^^,  xj^^,  xj^^  must  correspond  to  those  entered  for  these 

quantities  in  the  mesh  spacings  and  coordinates .  •  X^^^  cannot  be  greater  than 

F 

the  coordinate  corresponding  to  i  .  in  the  fine  mesh.  These  coordinates 

^  °  mm 

are  followed  by  values  for  the  relative  power  densities  (p^.^^.  , )  for  each  of 

i>  J»  k 

the  boxes  enclosed  by  the  coordinates  of  the  coarse  mesh  ..  The  values  of  the 


A  configuration  index  C  =0  indicates  that  the  entire  plane  is  of 
uniform  composition. 


relative  power  density  are  entered  in  the  following  order; 

„(0)  ^(0)  (0)  (0)  (0)  (0) 

^0.0,0’  ^0,0.1’  ••  ’  Po.O,K-l’  '^0,  1.0’  ■  '’Pq.I.K-I . ^0.J-1,K-1’ 

^1,0,0’  ■  ^i-i.j-i.K-r 

If  and  are  the  nunaber  of  coordinates  for  the  three  axes  then 

there  are  values  of  relative  power  density  listed. 
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IV.  GENERAL  PROCEDURE 

TRIXY  can  be  divided  into  six  sections 

1 .  Master  Tape  Generation 

2.  Rough  Source  Guess 

3.  Group  Source 

4.  Group  Sweep 

5.  Fission  Source 

6.  Edit  Code 

The  relation  between  the  various  sections  is  illustrated  by  the  overall 


TRIXY  flow  chart  (see  Figure  3) . 


OVERALL  FLOW  CHART 
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V.  GENERATION  OF  THE  MASTER  TAPE  AND 
THE  ROUGH  FISSION  SOURCE  GUESS 


This  section  of  TRIXY  consists  of  three  separate  programs . 

Master  tape  1  processes  and  checks  the  input,  generates  the  I  planes 
of  composition  indices,  and  generates  a  working  library  tape. 

Master  tape  2  generates  the  master  tape  for  all  groups  up  to  the  Fission 
Source  program. 

Master  tape  3  finishes  the  master  tape  and  generates  the  Rough  Fission 
Source  Guess,  if  needed. 


EQUATIONS  -  MASTER  TAPE  GENERATION 


Backward  Coefficients  - 

Current 

Boundaries 

M)8  - 

0,j.k 

z(i)g 

1 

21 
m,  n=0 

A<2) 

j-m 

A<3) 

k-n 

/4 

(V  1) 

r(2)g  _ 

i.0,k  ~ 

z(2)g 

1 

^i-i 

k-n 

> 

/4 

(V  2) 

i,  n=0 

r(3)g  _ 

i.j.O 

z(3)g 

1 

i-l 

A<^> 

j-m 

U 

(V  3) 

5,m=0 

Forward  Coefficients  -  Current  Boundaries 


cV>f,  =  Z»>S  ^  aP  /4 

1,  j,k  I  ^ - -  j-m  k-n 


(V  4) 


21 


,(2)g 

i.J.k 

_  7(2)g 

1 

21 

L  n=0 

A<'> 

1-1 

(V  5) 

,(3)g 

1 

A<2>  /4 

(V  6) 

i.  j.  K 

K 

i,  m=0 

1-1 

j-m  / 

Backward  Coefficients  -  Flux  Boundaries  Substitute  i,  j,  or  k=l  in 


the  appropriate  formula. 


p(l)g  _ 
i.j.k 

1 

m,  n=0 

D®  .  ^ 

1,  j-m,  k-n 

A  (2)  A<3> 

j-m  k-n 

A 

r(2)g  . 

i.j.k 

1 

21 
£,  n=0 

D®  .  u 

i-A,].k-n 

h 

t-1  k-n  / 

A<^> 

J 

f.(3)g  _ 

i.j.k 

1 

d8 

i-C,  j-m,  k 

A<‘>  A<2)  / 

1-^  J-m  j 

4  (3) 

4  A, 
k 

i,  m=0 

Additional 

Coefficients 

V 


i,  j.  k 


Jl)g 

i,  j.k 


V  D®  . 


J2)g 
i,  j.k 


Z 


1=0 


D 


g 

i-1,  j.k 


A 


(1) 

i-fi 


(V  7) 


(V  8) 


(V9) 


(V  10) 

(V  11) 


1 


(V  12) 


1 


(V13) 


jg  _  ^g  y 

»  - h  i-i.j-m.k-n  i-H.j-m,  k-n 

I,  m,  n=U 


1 


- 

i>  j>  ^  0^ 


5,  m,  n=0 


i,  til;  n=0 


i,  j.k 


I,  m, n=0 

I.  J.  K 


M®,  .  ,  V  ,  , 

1-1,  j-m,k-n  i-i,  j-m,k-n 


F®  V 

i-l,j-m,k-n  i-l,k-m,k-n 


V  li 

k-n  i-1,  j-m,k-n 


=  I 


i,  j,k=l 


i,  j.k 


EQUATIONS  -  ROUGH  FISSION  SOURCE  GUESS 
Rough  Power  Guess 

1,  J.k 

Unnormalized  Fission  Source 

=  P®’,  tV.  . 

i,j,k  i,j,k  i,j,k 

Normalized  Fission  Source 


i.J.k  i,j,k' 


J.  K 


^(OH 

j,k=l  i,  j, 


(V  14) 

(V  15) 

(V16) 

(V17) 


(V18) 


(V  19) 
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MASTER  TAPE  1  (see  Figure  4) 

Input  IS  read  into  the  machine.  The  program  is  removed  from  the 
machine  and  the  input  corrected  if  any  of  the  following  errors  occur; 

1.  The  Mesh  Spacings  and  Coordinates  do  not  agree  with  I,  J,  or  K. 

2 .  The  value  of  G  given  m  the  input  differs  from  that  in  the  file  selected 
from  the  cross  section  tape. 

3.  The  product  IxJxKxG  is  greater  than  150, 000. 

4 .  List  A  or  B  is  out  of  order . 

5.  The  configuration  is  outside  the  I  or  J  boundary. 

6.  The  product  J  x  K  is  greater  than  1000  for  8000  word  memory,  2000 
for  16, 000  word  memory,  or  4000  for  32, 000  word  memory. 

7.  Elastic  Scattering  was  requested  (this  is  not  coded  at  present). 

8.  There  are  too  many  words  of  input  in  the  Geometrical  and  Physical 
Data,  List  A,  List  B,  Mesh  Spacings  and  Coordinates,  Edit  Input,  or 
the  Rou^  Power  Guess . 

The  maximum  number  of  words  allowed  in  the  various  sections  of 
input  is  as  follows: 

Input  Maximum  No.  of  Words 

Geometrical  and  Physical  Data  50 

List  A 


List  B 


1500 

1000 
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Input  Maximum  No.  of  Words 

Mesh  spacings  and  Coordinates  250 

Edit  Input  550 

Roug^i  Power  Guess  500 

After  the  input  has  been  checked,  this  section  generates  k-planes  of 
composition  indices  and  writes  them  on  tape  three.  Each  i-column  is  written 
as  a  record  and  each  k-plane  as  a  file.  After  all  k-planes  have  been  generated, 
the  i -planes  are  generated  by  sorting  the  k-planes  and  are  written  one  record 
per  plane  on  tape  2.  It  is  possible  to  remove  the  program  from  the  machine  at 
this  time  so  that  control  rod  positioning  may  be  checked. 

The  program  constructs  a  working  library  of  the  remaining  programs 
on  tape  5  and  writes  the  Group  Source  program  on  tape  3.  The  cross  section 
records  are  read  from  tape  4,  and  Master  Tape  2  of  the  Master  Tape  code 

is  read  from  Tape  1.  Tapes  1  and  4  are  then  replaced  by  blank  tapes. 
MASTER  TAPE  2  (see  Figure  5) 

This  section  writes  part  one  of  the  Group  Sweep  restart  procedure  on 
the  master  tape  (or  tapes) .  This  record  is  followed  by  the  file  containing  the 
Group  Source  program . 

Three  passes  over  all  planes  in  the  mesh  are  made,  unless  the  plane  is 

(2) 

a  flux  boundary.  The  first  pass  computes  the  sentinals  and  C  backward 
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coefficients  which  are  written  on  the  master  tape  as  one  record. 

The  second  pass  computes  the  E,  r^^^,  coefficients .  These 

coefficients  are  written  as  records,  which  vary  in  length  depending  on  memory 
size  but  are  always  an  integral  number  of  columns .  Each  column  is  preceded 

by  a  block  of  r^^^,  r^^\  backward  coefficients .  The  coefficients  E,  r^^^ 

(2)  (3)  (1)  (2)  (31 

r'  C'  \  and  backward  coefficients  r  ,  r  ,  C'  '  are  omitted  if  the  point 

lies  on  a  flux  boundary  If  the  last  i-plane  is  a  current  boundary,  the  r^^^. 

A,  J>  X 

are  replaced  by  ,  .  Also,  if  k=K  is  a  current  boundary,  then 

I.J.k  i,j,K 

are  computed  by  the  formula  given  for  boundary  points .  If  the  plane  j=J  is  a 

(2)  (2) 

current  boundary,  r  are  replaced  by  C  . 

i>  J  >  k  i>  X 

(g)  ’(gi  "  (g) 

The  third  pass  is  made  to  compute  H  ,  F  ®  ,  and  V  .  H  are 

(g) 

computed  for  all  planes  unless  the  plane  is  a  zero  flux  boundary.  H  ®  are 

’(g) 

written  as  one  record  (JK  words)  per  plane  on  tape  1 .  The  F  are  written 

F  F 

on  tape  4  as  one  record  per  plane.  Only  the  planes  between  i  and  i 

F  F 

inclusive  are  written,  and  on  those  planes  only  the  points  between  j  .  ,  j 

min .  max , 

F  F  ’(g) 

and  k  ,  k  inclusive  are  written.  F  are  computed  only  for  groups 

min.  max. 

g  >  g  .  V  are  computed  only  during  the  mesh  pass  when  g=G.  are 

not  computed  for  this  group . 

At  the  end  of  a  group  g  <  G,  the  Group  Sweep  code  is  transferred  from 
tape  5  to  the  master  tape  following  a  file  of  E,  r^^^,  r^^\  The  Group 
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Source  code  is  written  after  the  file  containing  the  Group  Sweep  code.  The 
(sc) 

H  are  transferred  to  the  master  tape  from  tape  1 .  The  mesh  passes  are 
then  repeated  for  the  next  higher  group .  At  the  end  of  group  g=G,  the  Group 
Sweep  code  is  transferred  to  the  master  tape  and  the  code  for  Master  Tape  3 
is  read  into  memory. 

MASTER  TAPE  3  (see  Figure  6) 

This  section  transfers  the  Fission  Source  code  from  tape  5  to  the  (Figure  7 ) 

'(g) 

master  tape.  The  following  files  are  written  on  the  master  tape;  (1)  the  F 
one  group  per  file,  (2)  the  Edit  code,  (3)  Edit  input,  (4)  V  ,  (5)  mesh  spacings, 
and  (6)  the  I-planes  of  composition  indices .  The  Source  Flux  Guess  Word 
(SFGW)  is  inspected  to  see  if  a  flux  guess  tape  is  supplied  from  a  previous 
problem .  If  it  is  supplied,  the  machine  stops  and  the  flux  guess  tape  is  put 
on  tape  5.  If  a  rough  power  guess  is  supplied,  it  is  expanded.  The  unnormal¬ 
ized  source  guess  is  computed  and  written  on  tape  4.  This  guess  is  then 
normalized  and  written  on  tape  5 .  The  block  of  input  needed  by  the  Group 
Source  and  Group  Sweep  programs  is  set  up.  The  Group  Source  code  is  read 
in  from  tape  3 . 

It  is  now  possible  to  provide  a  duplicate  of  the  master  tape  if  one  was 
not  written  simultaneously  with  the  first  cc^y.  After  duplication  of  the  master 
tape  (if  performed),  control  is  transferred  to  the  Group  Source  program. 
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TAPE  UNIT  OPERATION 

Every  record  that  is  written  on  tape  is  backspaced  and  read  to  check 
for  errors.  If  there  is  a  tape  check  or  a  check  sum  error,  the  record  is  re¬ 
written  and  reread.  If  there  is  still  an  error,  the  machine  stops ,  When 
reading  a  record,  the  check  sum  and  tape  check  are  tested;  if  either  is  in 
error  the  record  is  backspaced  and  reread.  If  there  is  still  an  error,  the 
machine  stops .  If  a  tape  error  occurs  in  master  tape  1 ,  the  program  must  be 
reloaded  and  started  again.  If  the  error  occurs  in  master  tape  2  and  master 
tape  3,  it  is  possible  to  try  again  to  rewrite  or  reread  the  record.  If  the  error 
repeats,  the  entire  program  must  be  restarted. 

The  records  on  tapes  2,  3,  4,  and  5,  are  numbered.  The  program 
checks  to  see  that  it  has  the  correct  record.  If  the  tape  is  mispositioned,  the 
program  backspaces  or  skips  over  the  proper  number  of  records  and  checks 
again  to  see  if  it  has  the  correct  record.  If  the  correct  record  is  not  found, 
the  machine  stops.  In  master  tape  2  and  master  tape  3  it  is  possible  to  try  to 
read  the  record  once  more .  If  the  correct  record  cannot  be  found,  the  entire 


program  must  be  restarted. 


MASTER  TAPE  PROGRAM  -  1 
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Figure  4 
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Figure  5  (2) 
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Rough  Fission  Source  Guess 
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Figure  7 
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VI.  GROUP  SOURCE  PROGRAM 

PURPOSE 

The  purpose  of  the  TRIXY  Group  Source  Program  (GSP)  is  to  prepare 
certain  data  for  use  by  the  Group  Sweep  Program.  These  data  consist  of  a  file 
of  SErrC  and  associated  coefficients  and  a  file  of  <p  (flux)  values  over  which 
iteration  is  to  be  performed.  In  addition,  the  Group  Source  Program  saves  on 
tape  the  results  of  all  previous  flux  iterations  so  that  they  are  readily  available 
either  to  the  Editing  Program  (in  case  of  convergence)  or  for  a  subsequent  run 
(in  case  of  non -convergence).  (See  Figures  8  and  9.  ) 

INPUT 

To  allow  these  purposes  to  be  accomplished,  the  following  information 
is  written  on  tape  by  previous  programs :  (See  Figures  10,  11,  and  12), 

1,  Tape  A: 

a.  Files  of  ErrC,  associated  coefficients,  and  sentinels  for  each 
group .  The  history  (identification)  record  of  each  ErrC  file 
contains  value  of  JC  (fraction  of  the  fission  spectrum  occurring 
in  each  group)  used  by  the  Group  Source  Program  in  calculating 
S  (group  source).  It  also  contains  other  values  which  are  trans¬ 
ferred  to  the  history  records  of  the  files  for  use  by  the 


Group  Sweep  Program . 


36 


b.  Files  of  H  (scattering  coefficients)  for  each  group  except  the 
last. 

2.  je  B: 

a .  A  file  of  iterated  for  the  previous  group  (not  applicable  when 

the  Group  Source  Program  is  in  group  1). 

3.  Tape  E: 

a.  A  file  of  8  (fission  source). 

b.  A  file  or  files  of  4>  guess.  (These  files  are  not  present  at  the 
start  of  the  first  run  if  a  ^  guess  is  made  by  the  Group  Source 
Program) . 

NOTE;  There  are  files  of  information  present  on  the  tapes  other  than 
those  enumerated  above.  However,  these  files  are  not  used  by  the 
Group  Source  Program . 

OPERATION 

The  operation  of  the  Group  Source  Program  can  be  described  by  the 
following  chronological  steps .  It  is  assumed  that  each  file  written  on  tape 
includes  a  history  record. 

1.  (Not  applicable  to  group  1).  A  plane  of  H  for  the  previous  group  is 
read  from  tape  A.  The  corresponding  plane  of  iterated  sp  is  read 
from  tape  B.  This  latter  plane  is  transferred  to  tape  E  (for  ultimate 
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use  by  the  Fission  Source  Program),  following  the  last  file  previously 
written  on  tape  E .  Points  of  H  are  multiplied  by  corresponding  points 
of  ^  ,  and  a  plane  of  H<^  is  written  on  tape  C .  This  entire  step  is 
repeated  for  all  planes,  and  an  end  of  file  is  written  on  tapes  C  and  E. 
Tapes  B,  C,  and  E  are  rewound. 

2.  The  values  of  J)C  pertinent  to  the  group  under  consideration  are  obtained 
from  the  ErrC  history  record  on  tape  A  and  stored  for  future  use. 

3.  A  plane  of  H  ^  is  read  from  tape  C.  If  a  record  of  backward  coefficients 
is  present,  it  is  read  from  tape  A.  A  record  of  sentinels  is  read  from 
tape  A  and  transferred  to  tape  D.  These  sentinels  are  also  retained  in 
memory  for  future  use.  A  plane  of  S  is  read  from  tape  E  and  value  of 

S  at  each  point  in  the  plane  is  calculated  from  the  formula 
S  =  j:  «  +  Hfp 

Appropriate  short-cuts  are  taken  in  this  calculation  in  case  (a) 

jC=  0  or  1,  (b)  2  is  not  given,  as  for  points  outside  the  fissionable 

region,  or  (c)  H  ^  is  not  present,  as  in  the  case  of  group  1 . 

4 .  A  block  of  ErrC  coefficients  is  read  from  tape  A .  Under  control  of  the 
sentinels,  the  ErrC  coefficients  are  written  on  tape  D  so  that  each  value 
of  E  is  immediately  preceded  by  the  corresponding  value  of  S .  The 
records  thus  produced  are  known  as  SErrC  records ,  At  this  point,  if 
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a  ^  guess  for  the  current  group  is  not  present  on  tape  E,  the  program 
proceeds  to  step  5.  If  a  guess  is  present,  a  test  is  made  to  see  if 
the  current  SErrC  plane  is  completed.  If  it  is  not  completed,  step  4  is 
repeated.  If  the  plane  is  completed,  the  program  proceeds  to  step  6. 

5.  A  flux  guess  is  made  for  each  point  in  the  current  block  according  to 
the  formula 

S/E. 

A  test  is  made  to  see  if  the  plane  is  completed.  If  the  plane  is  not 
completed,  step  4  is  repeated;  otherwise,  the  program  proceeds  to 
step  6. 

6.  A  test  is  made  to  see  whether  all  planes  are  processed.  If  they  are  not 
processed,  the  procedure  is  repeated  from  step  3.  After  all  planes  are 
completed,  an  end  of  file  is  written  on  tapes  B  and  D,  and  on  tape  E,  if 
a  guess  was  made  by  the  Group  Source  Program.  Tapes  B,  C,  and 
D  are  rewound.  The  <f>  guess  tape  (tape  B  if  a  guess  was  made, 
tape  E  if  it  was  given)  is  positioned  to  the  beginning  of  the  ^  file  over 
which  the  Group  Sweep  Program  is  to  iterate .  The  latter  program  is 
read  into  memory  and  control  is  transferred  to  it . 

FINAL  RESULTS 


After  all  groups  are  processed  for  a  run,  the  information  written  on 
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tapes  B,  C,  and  D  by  the  Group  Source  Program  is  erased  by  either  the  Group 
Sweep  Program  or  the  Fission  Source  Program.  However,  the  pertinent  data 

is  written  on  tape  E  for  use  by  the  Editing  Program  in  either  of  two  orders . 

1.  First  run  with  (p  guess  not  given: 

a)  2 

b)  £r  guess  for  group  1  written  by  G.S.P.  during  group  1 

c)  Iterated  0  for  group  1  written  by  G.S.P.  during  group  2 

d)  4>  guess  for  group  2  written  by  G.S.P.  during  group  2 

e)  Iterated  0  for  group  2  written  by  G.S.P.  during  group  3 

f)  0  guess  for  group  3  written  by  G.S.P.  during  group  3 

g)  Iterated  0  for  group  3  written  subsequently  by  Fission 

Source  Program 

2.  First  run  with  0  guess  given,  and  all  subsequent  runs: 

a)  % 

b)  f  (Absent  during  first  run) 

c)  0  guess  for  group  1 

d)  0  guess  for  group  2 

e)  0  guess  for  group  3 

f)  Iterated  0  for  group  1  written  by  G .  S .  P .  during  group  2 

g)  Iterated  0  for  group  2  written  by  G.S.P.  during  group  3 

h)  Iterated  0  for  group  3  written  subsequently  by  Fission 

Source  Program 
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Figure  11 


HISTORY  RECORD  SUBROUTINE 


49 


50 


VII.  GROUP  SWEEP  PROGRAM  * 

The  TRDCY  Group  Sweep  Progpram  improves  the  fluxes  of  each  energy 

group  (g)  by  applying  the  "accelerated  Liebmann  technique.  "  The  flux 

estimates  are  given  at  mesh  points  i,  j,  k  and  are  improved  as  follows.  Each 

(3) 

point  along  a  line  in  the  X'  (or  K)  direction  (i,  j  constant)  is  improved  in 
successive  order  from  k=l  to  k=K;  this  process  is  called  a  "line  sweep . "  In 
a  given  plane  (i  constant),  the  line  sweeps  are  made  in  successive  order  from 
j=l  to  j=J,  and  this  process  is  referred  to  as  a  "plane  sweep."  Similarly,  a 
"mesh  sweep"  is  defined  as  the  process  of  making  plane  sweeps  in  successive 
order  from  i=l  to  i=I.  In  theory,  a  group  iteration  consists  of  as  many 
sweeps  as  are  necessary  to  satisfy  the  following  predetermined  condition; 

^  4>  ’  (vn  1) 

where  S  ^  is  the  primary  convergence  criterion  for  the  fluxes.  However,  two 
other  means  for  terminating  the  group  iteration  are  used: 

€^0  <  Cj%  of  (VII  2) 

where  is  the  percentage  of  initial  total  improvement  that  serves  as  a 

secondary  convergence  criterion,  and 


t  >  t 
-  max 


(VII  3) 


where  the  upper  limit  on  the  number  of  sweeps  allowable,  is  the 

*  Group  Sweep  Flow  Charts  (Figure  13) 
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g(t-l) 


tertiary  "convergence  criterion."  For  the  point  iteration  to  improve  .  i. 

ii  ^ 


g(t) 


the  following  fluxes  at  the  six  neighboring  points  are  used; 

^g<t)  ^g(t)  0  ,  and  .  After  ea, 

^i,j-l,k’  i,j,k-l’  i+l,j,k  i,  j+l,k  i,  j,k+l 

point  iteration,  the  improved  value  is  used  for  the  remainder  of  the  mesh 


sweep.  The  "Liebmann  Acceleration  Factor"  used  in  the  point  iteration 
equation  is  constant  throughout  a  group  iteration 

An  additional  factor  that  must  be  considered  is  the  neutron  balance  at 
the  end  of  each  mesh  sweep .  It  is  taken  into  account  by  a  renormalization 
factor,  Moreover,  the  flux  change  between  successive  iterations, 

g  is  calculated  at  each  mesh  point,  and  the  neutron  leakage,  L 

etc , ,  is  computed  for  each  mesh  point  on  a  boundary  plane . 

ASSUMPTIONS 


The  Group  Sweep  Program  is  stored  on  a  master  tape,  read  into 
memory,  and  given  control  by  the  Group  Source  Program .  The  following 
quantities  are  available  in  a  block  of  core  storage  (see  "Glossary  of  Symbols"); 
Location  (Decimal)  Decrement  Tag  Address 

29  at  (in  floating  point  form) 

30  g’  0  G 

31  g^  C  RCW 


32 


0 


r 
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Location  (Decimal) 

Decrement 

Tag 

Address 

33 

-- 

0 

J' 

34 

-- 

0 

K’ 

35 

-- 

0 

I 

36 

-- 

0 

J 

37 

-- 

0 

K 

38 

tape  A 

0 

tape  B 

39 

tape  C 

0 

tape  D 

40 

tape  E 

0 

flux  guess  tape 

41 

run  number 

0 

tape  F 

42 

0 

4 

1  (in'Jl,4) 

Tapes  are  supplied  as  follows: 

1 .  Tape  A:  master  tape  (positioned  to  read  Group  Source  Calculation 
Routine) . 

2 .  Tape  B:  available  for  working  storage  (rewound)  or  the  flux  guess  tape 
(see  5.  below). 

3.  Tape  C:  available  for  working  storage  (rewound), 

4.  Tape  D;  coefficient  tape  containing  the  following  file  (rewound): 
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Record  Number 
0 

1 

2 


3 


4 

N 


Contents 


History  record 

Plane  (J’K*  values)  of  ,  coefficients .  ^ 

K 

Word  count  (in  decrement),  J’  (or  2J')  sentinels; 
(2) 

K'  values  of  C  ,  ,  coefficients 

J-l,k 

Column  count  (in  decrement)  and  word  count 


(in  address); 


r<^> 

i-J.  j.k’ 


r(2)  AS) 

i.j-l.k'  i,j.k-l 


;  SB  r 


(1) 

i,  j,k' 


,  for  an  integral  number  of  columns 
i,j,k’  i,j,k 


Successive  records  of 
3  until  all  values  for 
plane  i  have  been 
recorded 


N+1  Record  number  2  for  plane  i+1 

N+2  Repeat  of  record  number  3  for  the  i+1  plane 

This  pattern  is  repeated  until  all  I'  planes  have  been  recorded. 
The  last  record  is  followed  by  an  end  of  file  gap , 

5.  Flux  guess  tape;  tape  E  or  B  depending  on  whether  or  not  a  flux  guess 
is  available  from  a  previous  problem,  positioned  to  the  file  containing 
the  following: 

^  This  record  is  not  present  if  plane  i=l  is  a  zero  flux  plane. 
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Record  Number 


Contents 


0 

1 

2 

N 


History  record 

Plane  (JK  values)  of  flux  estimates  for  plane  i=l^ 

Successive  records  of 
i -planes  until  all 
the  flux  guess  have 
been  recorded 


In  general,  N=r  .  However,  if  either  planes  i=l  or  i=I  are 
non  zero  flux  planes,  then  N=I. 

6.  Tape  F:  available  for  working  storage  (rewound). 

PROGRAMMING  METHOD 
Point  Iteration 

For  the  iteration  at  the  point  i,  j,  k,  three  coefficients  are  available 
from  previous  point  iterations .  (Boundary  points  are  exceptions  and  are 
discussed  later.)  The  coefficient  .  is  held  in  memory  from  an 

i”  K 

(2) 

iteration  in  the  previous  plane.  C'.  .  ,  v.  is  held  in  memory  from  an  itera- 

1>  J “  A? 

tion  in  the  previous  column  (a  "column"  is  a  line  of  mesh  points  along  the 

(3)  /3\ 

X'  '  direction).  C\  .  .  is  retained  from  the  previous  point  iteration. 

1,  J,  K-  i 

During  a  plane  sweep  of  the  group  iteration,  fluxes  for  the  planes 
i-1,  i,  and  i+l  are  held  in  memory;  this  requires  3JK  words.  Values  of 


This  record  is  not  present  if  plane  i=l  is  a  zero  flux  plane. 
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.  held  in  memory  require  JK  words,  while  retention  of 
requires  K  additional  words .  Sentinels  that  are  described  below  require 
J  words  if  K  <36  or  2J  words  if  K  >  36. 

Sources  and  Coefficients 

Sources  and  coefficients  are  read  into  their  reserved  memory  locations 
from  tape  D,  one  record  at  a  time.  Each  record  of  sources  and  coefficients 
contains  values  of  "SErrC"  that  have  been  calculated  prior  to  the  group  iter¬ 
ation.  For  each  mesh  point,  the  record  contains  either  the  set  of  two  values, 

S.  .  ,  and  E.  ,  or  the  set  of  five  values,  S.  .  , ,  E.  ,, 

i.j.k  i.J.k  i,j,k  i,j,k  i.j,k  i,j,k 

(3) 

and  C.  .  ,  .  The  sentinels  are  used  to  indicate  which  of  these  conditions  holds . 
i»  J*  K 

If  the  sentinel  bit  for  a  given  k  is  "O",  S  and  E  are  given;  if  the  bit  is  "1” 

the  set  of  five  values  is  supplied  .  The  sources  and  coefficients  are  stored  in 

mesh  sweep  order;  the  three  "backward”  coefficients,  ,  , , 

1-1, j.k’  i,j-l.k’ 

(3) 

and  C.  ,  are  stored  at  the  beginning  of  each  "column  '  of  SErrC  information , 

1,  J,  K- 1 

After  the  Group  Sweep  Routine  and  the  quantities  mentioned  in  the  preceding 
paragraph  have  been  stored,  the  memory  space  available  determines  the 
number  of  words  in  each  SErrC  record  for  a  particular  case.  If  possible, 
the  record  contains  the  sources  and  coefficients  for  an  entire  plane;  if  not, 
the  record  contains  SErrC  information  for  an  integral  number  of  columns ,  ^ 


Each  record  contains  a  word  count  and  a  column  count. 
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Boundary  Conditions 


Boundary  planes  i=l,  j=l,  and  k=l  are  exceptions  to  the  procedure 
described  above,  since  the  coefficients,  ,  and  . 

are  not  available  from  previous  points .  This  problem  is  resolved  by  calculat¬ 
ing  the  missing  coefficients  prior  to  the  group  iteration  and  storing  them  on 
tape  D.  The  first  record  after  the  history  record  on  tape  D  is  a  plane  of 

.  .  values;  however,  if  i=l  is  a  zero  current  boundary,  then  the  plane 

i“ii J*  K 

of  .  ,  values  is  set  to  zero  by  the  Group  Sweep  Routine .  There  is  also  a 

1"  ]» 

(2) 

column  of  C.  '  .  coefficients  and  the  sentinels  as  a  record  preceding  the 

(3) 

SErrC  information  for  each  plane.  The  C'  '  coefficient  is  supplied  with 

i,  J,  K- 1 

r^,^\  ,  ,  and  ,  ,  as  the  first  set  of  coefficients  for  each  column  in  the 
i-l,j,k  l,j.k-l 

SErrC  record (s)  for  a  plane. 

Iterations  are  skipped  for  mesh  points  in  a  specified  flux  boundary 

plane.  Thus,  if  plane  i=l  is  a  flux  plane,  flux  iteration  begins  on  plane 

i=2.  The  first  plane  of  flux  guesses  supplied  is  ^  ■  y  ~  •  i,» 

i"i>  Ij  ]» ■*^ 

the  plane  of  "backward"  coefficients  is  C^3^.  , ,  etc.  This  problem  of  logic 

i>  J»  -K 

is  handled  by  examining  the  boundary  conditions  (as  described  in  the  Reactor 
Control  Word)  at  the  beginning  of  a  given  group  sweep  and  presetting  several 


forks  throu^out  the  routine  (see  "Flow  Diagram"  for  a  detailed  description  of 


this  procedure). 
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Note  that  for  any  current  plane,  the  iteration  equation  contains  fluxes 
for  points  that  are  outside  the  mesh .  To  save  space,  the  memory  locations  for 
these  fluxes  are  not  carried  as  zeros .  Instead,  certain  instructions  used  in 
computing  the  flux  iteration  equation  are  changed  to  ”no  operation"  instructions 
and  then  are  restored  to  their  original  form . 

Use  of  Tapes 

For  the  first  mesh  sweep,  flux  estimates  are  read  from  tape  B  (or  tape 
E)  and  are  written  on  tape  C  one  plane  at  a  tune.  At  the  end  of  the  first  sweep, 
the  tape  designations  are  interchanged  as  follows;  the  flux  guess  tape  (tape  B 
or  tape  E)  is  designated  as  tape  F,  the  tape  saved  for  the  Repeat  and  Restart 
Procedures  (see  "Repeat  and  Restart  Procedures"),  tape  C  is  designated 
tape  B;  and  tape  F  is  designated  as  tape  C.  During  the  second  mesh  sweep, 
and  for  all  succeeding  sweeps,  fluxes  are  read  from  the  tape  designated  as  B. 
After  the  second  sweep,  the  tape  designations  are  interchanged  as  follows: 
tape  C  is  designated  as  tape  B;  tape  B  is  designated  as  tape  F;  and  the  original 
tape  B  is  designated  as  tape  C,  After  the  first  sweep,  the  flux  guess  tape  is 
available  for  restart  procedures .  After  the  second  sweep,  nothing  is  written 
on  the  flux  guess  tape.  Consequently,  if  tape  E  is  the  flux  guess  tape,  it  is 
positioned  at  the  end  of  the  file  of  flux  guess  planes  and  is  not  rewound  or  used 
for  writing.  After  the  third  sweep,  and  for  all  succeeding  sweeps,  the  tape 
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designations  are  interchanged  cyclically;  that  is,  the  tape  just  read  becomes 
the  restart  tape,  the  tape  just  written  becomes  the  tape  to  read,  and  the  tape 
previously  saved  for  restart  is  now  available  for  writing.  The  use  of  relative 
tape  designations  permits  this  procedure.  The  converged  flux  planes  always 
appear  on  the  tape  designated  as  tape  B.  The  flux  planes  are  written  on  the 
tape  designated  tape  C,  but  when  convergence  is  reached,  tape  B  and  tape  C 
designations  are  interchanged 
Procedure  After  Convergence 

The  group  number  (g^)  of  the  flux  estimates  to  be  improved  is  supplied 
in  the  history  records  of  tapes  D  and  flux  guess  tape  and  in  a  memory  location 
(see  "Programming  Conventions ") .  This  group  number  is  rewritten  on  the 
history  record  of  tape  B  to  indicate  the  group  number  of  the  converged  fluxes 
contained  on  the  tape .  After  convergence,  the  group  number  (gp)  stored  in 
memory  is  compared  with  G,  which  is  also  stored  in  memory.  If  gP  <  G, 
g^  is  increased  by  one  (so  that  gp  =  +  1).  and  the  Group  Source  Program 

is  read  into  memory.  If  g^  =  G,  the  run  number  is  increased  by  one,  and 
the  Fission  Source  Program  is  read  into  memory.  In  either  case,  tapes  B,  C, 

D,  and  F  are  rewound  and  control  is  transferred  to  the  beginning  of  the  program 
read  into  storage. 


59 


Storage  Allocations 

At  the  beginning  of  the  Group  Sweep  Program  the  following  list  of 
addresses  is  computed  from  control  data: 

Master  List  of  Initial  Addresses 


Location 

ML 

ML+1 


Address 

MS 

MS+J'  +  1  =  (MS)^ 
or 

MS  +  2J'  +  1  =  (MS), 


Description 

Sentinels 


(2) 

C  .  .  ,  .  column 
1,  j-l,k 


ML+2 

(MS)^ 

+  K’ 

+ 

1  =  (MS)^ 

•  1  plane 

ML+3 

(MS)^ 

+  JK 

+ 

1  =  (MS)^ 

BLOCK  1  of  flux 

ML+4 

(MS)3 

+  JK 

+ 

1  =  (MS)^ 

BLOCK  2  of  flux 

ML+5 

(MS)^ 

+  JK 

+ 

1  =  (MS)^ 

BLOCK  3  of  flux 

ML+6 

(MS)5 

+  JK 

= 

(MS)6 

SErrC 

Here  "MS"  is  the  address  of  the  first  word  of  memory  available  for 
storage.  (The  Group  Sweep  Routine,  associated  subroutines,  temporary 
storage,  and  control  data  are  stored  in  the  lower  part  of  memory  and  MS  is  the 
location  of  the  first  word  after  this  block.) 

The  maximum  number  of  SErrC  that  can  be  retained  in  memory  at  one 
time  by  the  Group  Sweep  Routine  can  be  computed  by  the  following  formula; 

NsErrC  =  "^Machine  -  (MS  +  5  +  ]•  +  K'  +  4,K),  (VII  4) 
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where  N,,  .  .  is  the  number  of  words  of  core  storage 

Machine  ® 

in  the  machine,  J*  =  J'  if  K  <  36,  or  J*  =  2J’ if  K  >  36. 
This  formula  is  used  in  the  master  tape  routine  to  determine  the  maximum 
number  of  ErrC  columns  that  can  be  written  in  one  record.  If  the  Group  Sweep 
Routine  is  reassembled,  the  absolute  location  for  MS  may  change.  In  this  case, 
the  constant  value  of  MS  used  in  the  master  tape  routine  also  must  be  changed. 

Since  storage  locations  are  assigned  within  the  program,  increased 
memory  capacity  is  utilized  without  changes  in  programming.  The  limits 
imposed  on  input  quantities  I,  J,  and  K  control  storage  allocations . 

Tape  Subroutines 

All  reading  or  writing  on  tapes  B,  C,  D,  E,  or  F  is  done  with  a  tape 


read- write  subroutine.  When  writing,  the  subroutine  computes  a  complemented 
logical  check  sum  and  stores  it  as  the  last  word  of  the  record,  but  the  record 
is  not  backspaced  after  writing  to  test  the  check  sum.  However,  when  reading, 
the  check  sum  is  tested  as  well  as  the  tape  check  indicator,  and  if  either  test 
fails  the  subroutine  returns  to  the  error  return.  If  an  end  of  file  is  encountered, 
return  is  made  to  the  error  return. 

The  file  identification  records  of  programs  on  the  master  tape  are  read 
by  means  of  the  Read  File  Identification  Subroutine.  This  subroutine  reads  the 


file  identification,  the  first  word  of  the  program  history  record,  and  compares 
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it  with  the  file  identification  given  in  the  calling  sequence.  If  the  file  identi¬ 
fication  checks,  the  subroutine  reads  the  next  record  (the  program)  into 
temporary  storage  to  check  the  check  sum  and  the  status  of  the  tape  check 
indicator.  If  these  tests  are  all  ri^t,  the  master  tape  is  backspaced  one 
record,  the  Read  Program  Record  Routine  (RPR)  is  read  into  0  through  21 
(decimal),  and  the  subroutine  transfers  control  to  the  RPR  subroutine  (see  below). 
If  the  file  identification  does  not  check,  the  subroutine  positions  the  master 
tape  to  the  beginning  of  the  next  file,  checking  to  be  certain  that  the  end  of  the 
master  tape  has  not  been  reached,  and  reads  the  next  file  identification .  If 
the  proper  file  is  not  found,  this  process  continues  to  the  end  of  the  master 
tape  which  is  marked  by  a  double  end  of  file .  The  tape  is  then  rewound  and 
the  search  is  repeated.  If  the  correct  file  cannot  be  found  after  the  second 
pass,  the  subroutine  prints  the  error  indication  instructions,  calls  for  the 
extra  master  tape,  and  stops .  The  program  is  prepared  to  repeat  the  entire 
search  when  the  start  button  is  pressed.  If  either  the  check  sum  or  tape  check 
test  fails,  the  subroutine  repeats  reading  the  record  once  before  giving  the 
error  indication  as  above .  This  procedure  for  reading  the  file  identification 
is  an  attempt  to  assure  that  the  proper  program  is  found  and  read  into  memory, 
thereby  avoiding  a  stop  in  the  operation  of  TRIXY. 

The  RPR  subroutine  is  a  tape  read  routine  that  reads  the  program  into 
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core  storage,  beginning  at  location  transfers  control  to  The 

routine  reads  the  record  twice  before  giving  the  error  stop  at  This 

routine  must  be  relocated  to  0  through  2l  (decimi  is  used. 

Entering  New  Data 

During  the  operation  of  the  Group  Sweep  Program  it  may  be  desirable 
to  enter  a  new  Liebmann  acceleration  factor  {ai).  At  the  beginning  of  a  mesh 
sweep,  it  is  possible  to  load  decimal  data  from  the  card  reader  by  putting  the 
deck  for  the  data  load  program  MDL9C  and  decimal  data  in  the  card  reader . 

If  the  card  reader  is  readied  and  sense  switch  4  depressed,  the  Group  Sweep 
Program  loads  the  cards .  MDL9C  is  read  into  upper  memory  by  a  special  load 
routine  (MBILD)* and  destroyed  later  by  reading  in  flux  planes .  If  an  error  is 
found  in  the  data,  a  comment  is  printed  and  the  program  stops  at  (42)^^.  If  an 
end  of  file  card  is  encountered,  the  program  stops  at  (42)^^.  If  the  normal  end 
of  record  punch  is  given,  return  is  made  to  the  Group  Sweep  Program.  The 
Group  Sweep  Program  prints  instructions  and  stops  at  (^2)^^  to  allow  these 
instructions  to  be  followed. 

Itermediate  Printing 

The  Group  Sweep  Program  uses  the  TRIXY  Print  Subroutine  (PR)  to 

print  certain  information  during  a  group  iteration. 

*MBILD  is  loaded  into  locations  0-24  and  does  not  destroy  the  control 
words  stored  in  29  -  42. 
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At  the  beginning  of  the  first  mesh  sweep  the  following  information  is  printed; 

(1)  Run  and  Group  Number 

(2)  Convergence  Code  and  Sweep  Number 

(3)  Liebmann  Acceleration  Factor  (  ) 

^  0  (primary  convergence  criterion) 

(5)  Cj^  (secondary  convergence  criterion) 

(6)  ’'max  convergence  criterion) 

At  the  end  of  each  mesh  sweep,  the  sweep  number  (t)  is  compared  with  a  print 
sweep  counter .  The  print  sweep  regulator  (C2)  is  stored  in  this  counter  on  the 
first  sweep,  and  C2  is  added  to  the  counter  each  time  it  is  equal  to  the  sweep 
number .  If  the  sweep  number  is  greater  than  or  equal  to  the  print  sweep 
counter  and  if  sense  switch  3  is  depressed,  the  following  quantities  are  printed: 

(1)  Run  and  Group  numbers 

(2)  Convergence  Code  and  Sweep  Numbers 

(3)  Liebmann  Acceleration  factor  (  <jC  ) 

(4)  Flux  Change 

(5)  Error  ratio  (  ^  ) 

(6)  Cj%of 

When  convergence  is  reached  this  information  is  printed  regardless  of  the 


setting  of  sense  switch  3. 
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Repeat  Procedure 

Since  the  three  tapes  used  for  reading  and  writing  flux  planes  are 
interchanged  cyclically  between  mesh  sweeps  (see  "Use  of  Tapes"),  during 
any  sweep  (t)  the  flux  planes  from  the  previous  sweep  (t-1)  are  available. 
However,  there  are  no  flux  planes  available  from  the  previous  sweep  when  t=l, 
and  the  repeat  procedure  described  below  is  not  applicable  during  the  first 
sweep.  Also,  an  exception  must  be  made  when  t=2.  If  the  flux  guess  is  supplied 
on  tape  E,  during  the  second  sweep  tape  E  is  the  tape  theoretically  available  for 
restart.  Tape  E  would  have  to  be  rewound  to  use  it  again  but  the  Group  Source 
Program  assumes  that  tape  E  is  not  rewound.  The  special  procedure  followed 
when  t=l  or  t=2  (and  the  flux  guess  tape  is  tape  E)  is  described  later. 

If  an  error  is  encountered  in  reading  flux  planes  during  a  sweep,  the 
record  is  reread.  If  the  error  occurs  again,  the  program  transfers  to  the 
Repeat  Procedure.  The  Repeat  Procedure  tests  the  sweep  counter .  K  tj^l  and 
t;^2  (or  if  when  t=2  the  flux  guess  tape  was  tape  B),  the  Repeat  Procedure  prints 
that  repeat  and  restart  procedures  are  available,  prints  instructions,  and  stops. 
When  the  start  button  is  pressed,  the  routine  transfers  to  Prepare  Restart 
Procedure  (PRP) .  If  sense  switch  5  is  not  down,  PRP  returns  control  to  the 
Repeat  Procedure.  At  this  point,  one  is  subtracted  from  the  sweep  counter, 

R^^  is  stored  in  place  of  R^^\  the  designations  of  tape  F  (which  contains  the 
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flux  planes  for  sweep  t-1)  and  tape  B  are  interchanged,  tapes  B,  C,  D,  and  F 
are  rewound,  and  control  is  transferred  to  the  beginning  of  the  mesh  sweep .  In 
other  words,  if  during  sweep  t’  an  error  is  encountered  in  reading  from  tape  B, 
the  Repeat  Procedure  begins  the  iteration  on  sweep  t'  -  1.  Since  all  conditions 
for  the  iteration  on  sweep  t’  -  1  have  been  restored,  using  the  Repeat  Procedure 
produces  no  change  in  the  results  of  a  mesh  sweep . 

NOTE;  If  the  error  in  reading  flux  planes  on  sweep  t'  occurs  a  second 
time,  it  would  be  desirable  to  remove  TRDCY  from  the  machine.  (Some  mal¬ 
function  of  the  machine  with  regard  to  tape  reading  or  writing  is  indicated,  and 
this  error  prevents  TRDCY  from  continuing).  In  this  case,  when  the  stop  occurs 
which  signifies  that  a  repeat  procedure  is  necessary,  sense  switch  5  may  be 
depressed  before  the  s>art  button  is  pressed.  As  described  above,  a  transfer 
to  Prepare  Restart  Procedure  (PRP)  is  executed  (see  below).  The  Repeat 
Procedure  is  designed  to  compensate  for  a  single  isolated  tape  error,  while 
the  Restart  Procedure  is  used  in  event  of  a  persistent  error. 

If  t=l  or  t=2  and  the  flux  guess  tape  is  tape  E,  the  Repeat  Procedure 
prints  that  repeat  and  restart  procedures  are  not  available  and  stops  at  the 
special  location  to  indicate  error . 

Prepare  Restart  Procedure 


In  addition  to  the  Repeat  Procedure  described  abo^’^p  a  method  is 
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provided  for  interrupting  TRIXY  in  the  middle  of  the  Group  Sweep  Program, 
removing  TRDCY  from  the  machine,  and  restarting  the  program  after  the 
interruption.  The  Restart  Procedure  is  designed  for  either  voluntary  inter¬ 
ruption,  (e.g. ,  when  the  machine  is  requested  for  the  operation  of  other  pro¬ 
grams)  or  mandatory  interruption  (in  the  event  of  a  persistent  machine  error 
that  demands  correction).  At  the  beginning  and  the  end  of  a  mesh  sweep,  the 
Group  Sweep  Program  transfers  to  the  PRP  subroutine.  If  sense  switch  5  is 
up,  PRP  immediately  returns  control  to  the  program.  If  sense  switch  5  is  down, 
the  routine  checks  the  sweep  counter .  If  t=l  or  t=2  and  the  flux  guess  tape  is 
tape  E,  the  routine  returns  control  to  the  main  program. 

If  t^l  and  tj^2  (or  if  t=2  and  the  flux  guess  tape  is  tape  B),  PRP 
prints  instructions  and  prepares  the  restart  procedure .  The  designations  of 
tapes  B  and  F  are  interchanged,  and  the  flux  guess  tape  is  designated  as  tape  B. 
The  subroutine  tests  to  see  if  tape  A  is  designated  as  logical  tape  1  or  3 .  If 
tape  A  is  logical  tape  1,  PRP  does  not  modify  the  tape  designations,  as  tape  A 
may  be  loaded  using  the  load  tape  button.  If  tape  A  is  not  logical  tape  1  but 
logical  tape  3,  PRP  sets  tape  A  to  logical  tape  1 .  The  subroutine  finds  the 
tape  formerly  designated  as  logical  tape  1  and  sets  this  tape  to  logical  tape  3, 
thus  preventing  two  tapes  from  being  designated  as  logical  tape  1 .  PRP  punches 
out  the  contents  of  locations  30  through  42  and  stops  at  a  special  location  to 
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indicate  that  TRIXY  may  be  removed  from  the  machine . 

To  restart  TRIXY  only  the  following  need  be  done- 

(1)  The  tapes  must  be  replaced  exactly  as  they  were.  If  this  is  the  first 
time  TRKY  has  been  restarted,  the  master  tape  (formerly  logical 
tape  3)  must  be  set  to  tape  1,  and  the  tape  that  was  1  must  be  set  to 
tape  3. 

(2)  The  card  punched  out  by  PRP,  followed  by  the  deck  for  the  data  load 

routine  (MDL9C)  and  a  data  card  giving  the  proper  Liebmann  factor 
(  ),  must  be  placed  in  the  card  reader. 

(3)  Sense  switch  4  must  be  depressed  and  the  other  sense  switches  set  as 
before . 

(4)  The  load  tape  button  must  be  pressed. 

(More  detailed  instructions  for  restarting  TRIXY  are  given  in  the  "TRIXY 
Operating  Manual  Appendix  E.  ) 

When  the  load  tape  button  is  pressed,  the  self-load  loop  on  the  master 
tape  is  read  into  memory .  The  data  contained  on  the  card  punched  out  by  PRP 
is  read  into  memory  locations  30  through  42,  thereby  re  storing  the  quantities 
that  various  TRIXY  routines  assume  are  in  memory.  The  Group  Source  Pro¬ 
gram  is  read  into  memory,  and  a  sense  light  is  turned  on  to  indicate  to  the 
Group  Source  Program  that  is  is  a  restart.  Control  is  transferred  to  the  Group 
Source  Program  which  places  tape  E  in  the  required  position .  (K  a  flux  guess 
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was  available  from  a  previous  problem,  tape  E  will  be  positioned  as  if  the 
flux  guess  planes  had  been  read  by  the  Group  Sweep  Program.)  The  Group 
Source  Program  reads  the  Group  Sweep  Program  into  memory  and  transfers 
control  to  the  Group  Sweep  Program.  At  this  point,  all  the  conditions  of  core 
storage,  tape  positions,  and  control  at  the  time  TRIXY  was  interrupted  have 
been  duplicated,  except  that  the  Liebmann  factor  (  at)  is  not  available.  Con¬ 
sequently,  at  must  be  read  into  memory  by  the  Group  Sweep  Program .  (See 
instructions  in  (2)  and  (3)  above.) 

If  TRIXY  was  interrupted  in  the  middle  of  sweep  t’,  the  Restart 
Procedure  starts  the  iteration  on  sweep  t’  -  1,  There  may  be  a  slight  difference 
in  the  results  on  sweep  t'  after  TRIXY  has  been  restarted,  because  the 
renormalization  factor  (R  (t)  )  is  set  equal  to  one  by  the  Group  Sweep  Routine. 
However,  the  effect  of  this  change  is  negligible,  and  all  tests  of  the  restart 
procedure  indicate  that  virtually  no  change  in  final  results  is  noticable  after 


a  restart. 


EQUATIONS 


Point  Iteration  Equation; 


-P 


.(t) 


(1 

■’] 

i,  j.k 

cC 

s 

+ 

- 

®i.  j.k 

i.j.k 

1-1.  j.k 

i-1.  j.k 

<p(t) 

1,  J.k- 

1  i,  j. 

,k-l 

\  i.J.k 

1.  J.k 

i,  j+1. 

k  1.  J.k 

1.  J.k+w 

(VII  5 ) 


where  .  .  is  the  improved  flux  calculated  from  the  flux  estimate 

if  J>  it 

,  and  A  is  the  Liebmann  Acceleration  Factor. 

1,  J,k 

Renormalization  factor  ):^ 


R 


(t-1)  _ 


I,  J,  K 


=  Z 


i,  j,k 


,-i 

S  ^  2 


(VII  6) 

X 

(t-l) 


i=l,j=l,k=l 
where  the  neutron  leakage  from  the  mesh  is 

=  l‘/>  +  L»)  ^  Lf)  ^  L<^>  ^  Lf>  ^  4^) 


>  E.  .  ,  f>',‘  , 

‘'''y 


(VII  7) 


For  specified  current  condition 

].  K 


-  5 


i=l,k=l 


c(i)  ^(t-1)  ^ 
0,j,k  l,j,k 


(VTI  8) 


(1)  - 


J.  K 

z: 

j=l.k=l 


k  -^k 

1>  j>  k  I,  j,  k 


(VII  9) 
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3. 


4. 


5. 


For  specified  flux  condition 

=  ■''5^  c</>. 

l.j.k  V  2,j,k 


?> 


j=l,k=l 

(1)  M)  A(t-l) 

''  *jinc=i  V  '-‘.J’*' 


(t-l)  A 

i.j.ky’ 

^I.j.k 


etc. 


Neutron  conductance  coefficients  (C): 

a) 


c<'>.  ,  =  r»)  ,  ^  r(.».  ,  , 

i.j.k  i,j,k  i.j,k-l 


b)  C<'>  + 

i,j,k  i.j.k  i,j.k-l 

Neutron  depletion  coefficient  (B): 

B  -  E  +  +  C' 

i»  j»  k  i,  j,k  i>j>k  i,  j,  k  i^j,  k  i“l,j,k 


,(1) 


+  c<'>  ,  ,  +  c<3>.  , 

i,j-l.k  i,j,k-l 


Flux  Change  (  € 

4'’  =  ^  ( 

^  i=lj=l,k=l^ 


i,j.k  ij,k 


Error  Ratio: 


where  C 


(t-C^) 

V 

is  the  frequency  print  regulator . 


(Vn  10) 

(VII  1 1 ) 

(VII  12) 
(VII  13) 

(VII  14) 

(VII  15) 


(VII  16) 
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7. 


Convergence  criterion: 


a) 


Primary 

<P 


(g) 


(VII  1) 


Q 

where  S  is  predetermined  convergence  criterion . 


b) 


Secondary 

.(t) 


<  of 


{VII  2) 


c) 


where  is  predetermined  percentage  (expressed  as  a  decimal) 
Tertiary 


t  >  t 

max 


(VII  3) 


where  t  is  the  upper  limit  on  the  number  of  allowable  sweeps , 
max 


Svsf££^ 


Gkoop  SwetLP 


Figure  1  3  (4 ) 


Figure  13(5) 


J^Af^T/AL  L.SAKAC&  QaL-CULAT^ON 


igure 
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VIII.  FISSION  SOURCE  PROGRAM  " 


DESCRIPTION  OF  PROGRAM 


The  TRIXY  Fission  Source  Program  computes  a  new  fission  source  by 
an  iterative  procedure  If  the  fission  source  has  not  converged,  tape  E  contain¬ 
ing  the  iterated  fission  source  and  the  flux  iterations  from  the  last  run  are 
prepared  and  turned  over  to  the  Group  Source  program,  unless  it  is  decided  to 
stop  the  fission  source  iterations. 

After  tape  E  has  been  prepared,  the  program  prints  out  several  lines  of 
information  which  include  the  criticality,  maximum  delta  S  and  its  location 
within  the  mesh,  the  total  fission  source  error,  and  the  value  of  used  for  the 
next  iteration.  Convergence  or  non- convergence  is  indicated.  The  program 
stops,  and  the  operator  can  go  directly  to  the  Edit  program  if  instructed, 

The  Fission  Source  Program  copies  the  iterated  fluxes  for  group  G  onto 


the  tape  containing  the  iterated  fluxes  for  the  previous  groups  and  the  fission 

'jT 

source  iteration  from  the  last  run.  Using  F.  .  .  from  ihe  master  tape  and  the 

g 

iterated  .  from  the  last  run,  the  following  values  are  computed  and 

1.  j,  Jc 

written  on  a  working  tape 


g*’’*  -  ^  or  CT 

j  6  _  p  o  o 

*i,j,k  i,j,k  'i,j,k 


(VIII  1) 


S  ®  ,  for  the  next  higher  group  are  computed,  and  i  ^  from  the 

i,j,k  1,  j,k 


-Fission  Source  Flow  Charts  (  Figure  15) 
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previous  run  are  read  into  the  memory  and  added  to  these  results.  The  sum  is 


* 

written  on  tape  and  the  process  repeated  until  the  sum  ^  ^  been  formed. 

For  each  group,  the  program  calculates 

F 


fS  =  *5"  X  /V 


(Vlll  2) 


and  the  sum 


f  =  ^  f^. 
g=g' 


When  all  groups  are  finished,  the  normalized  fission  source  j  k  ) 


(Vlll  3) 


j(r-l) 


for  the  last  run  and  the  unnormalized  fission  source  (  %  ,  )for  the  present 

\  i.j.k; 

run  are  read  into  memory.  The  normalized  source  is  computed: 


o  (r)  (r- 1) 

i.J.k 


/  o  G  • 

^  ^.j,k 

f 


.(r-1) 
^i,  j.k 


(VIII  4) 


The  value  for  <d  is  read  in  from  the  card  reader.  As  the  new  X  Is  computed, 
the  maximum  delta  ^  is  found,  the  total  fission  source  error  is  accumulated, 
and  the  convergence  is  checked.  The  new  S  is  written  on  the  new  tape  E,  and 
the  iterated  fluxes  from  the  previous  run  are  transferred  to  this  tape.  The 
program  prints  the  information  mentioned  earlier  and  stops.  (See  Figure  14). 


&8 


END  OF  RUN  1  NOT  CONVERGED 

TAPE  2  IS  THE  NEW  TAPE  E. 

DEPRESS  SWITCH  6  IF  EDIT  IS  DESIRED 
PUSH  START  BUTTON  TO  CONTINUE 


CRITICALITY  1.1186704 

TOTAL  SOURCE  ERROR  1121.6346 

MAX  DELTA  S  DURING  RUN  134.94965 
PLANE  3.  COLUMN  2.  POINT  2. 

OMEGA  1.0000000 


Example  of  Printing  by 
FISSION  SOURCE  PROGRAM 


Figure  14 
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Figure  15  ( 1 ) 
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IX.  EDIT  PROGRAM 


PURPOSE 


The  Edit  Program  selects,  sorts,  and,  where  necessary,  computes 
requested  output.  The  following  equations  are  used: 

Power  density  at  point  i,  j,  k 


if  j»  k  1,  j,  k  i,  j,  k 

Power  integral  for  column  i,  j 
.  F 


max 


1.3 


Z 


k=k^. 

min 


i.j.k  ' 


(IX  1) 


(IX  2) 


Power  density  for  column  i,  j 

F 


max 


P.  .  = 

i.J 


If] 


Z 


ijfk  * 


k=k 


min 


Flux  integral  for  composition  region  1 


n 


g 


n 


1,  ].  K 


i,  j,  k=l 


g 

i,  jf  k 


1 

Z 


(1X3) 


V  .  8 

I'-*  .J-m.k-n  p.  ,  o,  ^  , 

.£,m,  n-0  n  i-5,  j-m,  k-n 

(IX  4) 


*Edit  Flow  Charts  {Figure  18) 
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where 


n  ^  i,  j,  k 


jg  ’  p 

n  "i.j.k 


1  if  i  =.« 

n  i,j,k 

(IX  5) 

0  if  i  . 

n  i,j,k 

(IX  6) 

Average  flux  for  composition  region  ^  ^ 


Z 


n/  i,  j,k=l  i,m,  n=0 


i-il ,  j-m,  k-n 


S 


Q  ,  H 

n  i-l,j-m,k-n 


(1X7) 


Output  is  divided  into  seven  categories.  (See  Figure  16). 

Data  1:  P.  .  .  each  point  of  column  i,  j,  in  order  of  increasing  k, 
i?  k 

followed  by  ^  .  and  P 

i.J  i,} 

Data  2:  2  _  for  each  column  of  plane  i,  in  order  of  increasing  j. 

Data  3:  ^  for  each  column  of  plane  i,  in  order  of  increasing  j , 

Data  4;  Number  of  fissions  per  source  neutron  (f^)  for  all  groups, 

followed  by  the  total  number  of  fission  per  source  neutron. 

Data  5;  ^  ®  for  group  g  for  each  composition,  in  order  of  increasing 

n 

index  number  ji 

n 

—  g 

Data  6:  for  group  g  for  each  composition,  in  order  of  increasing 


n 

index  number  ^ 


n 
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a 

Data  7:  Flux,  i.  for  group  g,  for  each  point  of  column  i,  j,  in 

ij  J» 

order  of  increasing  k. 

The  general  flow  of  the  Edit  Program  is  as  follows .  A  plane  of  X  .  .  . 

Ji  K 

•  t 

values  and  a  plane  of  V.  .  ,  values  are  read  from  tapes  E  and  A  respectively . 

K 

Data  1,  2,  and  3  are  found  for  a  column.  Data  1  are  written  on  tape  B,  Data  2 
on  tape  C,  and  Data  3  on  tape  D.  When  the  above  procedure  has  been  repeated 
for  all  columns  and  for  all  planes,  the  information  on  tapes  C  and  D  is  trans¬ 
ferred  to  tape  B,  Data  4  are  then  transferred  from  tape  E  to  tape  B.  Mesh 
spacings  and  material  indices  for  a  plane  are  read  from  tape  A,  and  the 
corresponding  values  of  flux  are  read  from  tape  E.  The  computation  of  ^  ® 

and  of  is  begun,  and  Data  7  are  found  for  one  column  at  a  time  and 

n 

written  on  tape  D.  When  Data  5  are  completed  they  are  written  on  tape  B,  and 
Data  6  are  written  on  tape  C.  Data  5,  6,  7  are  found  for  the  remaining  groups, 
and  the  information  on  tape  C  and  D  is  transferred  to  tape  B.  The  information 
on  tape  B  is  read,  converted  to  BCD,  and  written  on  tape  10 .  Unrequested 
data  is  not'Computed,  and  when  no  sorting  is  necessary  data  are  written  directly 
on  tape  B.  Each  section  of  output  has  a  two  line  heading  containing  the  data 
type  (q)  ,  the  plane  (i),  the  column  (j),  the  group  (g),  the  run  (r),  the 
problem  number  and  the  cross  section  data  file  number. 

The  output  routine  GL  OUT2  (SHARE  distribution  number  84)  is  used. 
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It  has  been  modified  to  set  up  a  72  character  line  and  to  use  the  SHARE  standard 
printer  board  number  2,  The  output  routine  writes  numbers  in  the  format 
-YY  XXXXXXXX.  The  decimal  point  is  assumed  to  be  to  the  left  of  the  man- 
tissa.  Positive  signs  are  indicated  by  blanks,  and  leading  zeros  in  the  charact¬ 
eristic  are  not  printed.  If  the  mantissa  is  zero,  no  characteristic  is  printed. 

RE  EDIT 

If  the  master  tapes  A  and  E  are  saved,  a  re-edit  for  information  not 
requested  during  the  edit  may  be  performed. 

Operating  instructions ; 

Load  and  ready  tape  A  on  unit  1 . 

Load  and  ready  tape  E  on  unit  5 . 

Load  and  ready  working  tapes  2,  3,  and  4. 

Set  up  output  tape  10 . 

Place  deck  in  read  hopper  containing 
RE-EDIT  DECK  (6  cards) 

Edit  Input  Request  cards 

Do  ^t  put  blank  cards  on  the  back  of  the  deck. 


CLEAR  -  LOAD  CARDS . 


Error  stops: 

(0463l)„  Illegal  double  punch  on  input  cards.  Correct  card  and 

O 

Start  over . 

(11646)  Tape  error  in  reading  history  of  edit  code.  Depress 

8 

START  to  try  again. 

(04563)  ERROR  IN  EDIT  CODE  printed  on-line.  Depress  START 

8 

to  try  again. 

EDIT  REQUEST 

There  are  seven  different  items  of  output  which  may  be  requested. 
These  items  are  combinations  of  the  following  data,  where  the  output  data 
are  given  in  mesh  sweep  order. 

Data  1 :  for  all  i  planes . 

Power  density  (P.  .  .  )  for  each  point,  followed  by  the  power 

if  jf  K 

integral  (  ^ .  .)  and  the  power  density  (P.  .)  of  the  column. 

1,  J  1, J 

Data  la:  for  selected  columns . 

Same  as  Data  1 . 


Data  2:  for  all  i  planes. 


Power  integral  (8  ^  j)  for  every  column. 


Data  3:  for  all  i  planes. 


Power  density  for  every  column. 
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Lata  4: 


Data  5: 


Number  of  fissions  per  source  neutron  (f^),  for  all  groups  followed 
by  the  total  number  of  fissions  per  source  neutron . 
for  all  groups . 

Flux  integral  (  ^  ^  )  for  each  composition,  in  order  of 


n 


increasing  l  . 

n 

Data  6:  for  all  groups  . 

Average  flux  for  each  composition  in  order  of  increasing  ! 

Data  7:  for  group  g  for  all  i  planes. 

Flux  (  <P  ®  .  . )  for  each  point. 
i> 

Data  7a:  for  group  g  for  selected  columns. 

Same  as  data  7. 

The  following  table  indicates  which  data  are  included  unde:  the  various 
items  of  output.  Improper  combinations  of  items  are  also  indicated. 


DATA  Excludes  Items 


Item 

A1 

A2 

A3 

A4 

A(44g) 
B(i.  j) 
Cg(i.l) 


1 

lA 

2 

3 

4 

5 

6 

7 

7A 

A1 

A2 

A3 

A4 

A(44g) 

B 

c 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Figure  16 


98 


Requests  for  items  A1  through  A(4-(-g)  are  submitted  in  the  form  of  an 
octal  word  WORD  .  WORD  I  is  first  represented  as  a  binary  number  by 
placing  a  one  in  the  box  under  each  requested  item  and  a  zero  in  the  box  under 
each  item  not  needed.  All  the  boxes  must  be  filled  regardless  of  the  number  of 
groups  in  the  problem  ,  Each  set  of  three  binary  digits  represents  one  octal  as 
shown  below. 


First  3  binary  digits  Subsequent  octal  digits 


000  =  -10 

0 

001  =  +l 

1 

010  =  +2 

2 

on  =  +3 

3 

100  =  -0 

4 

101  =  -1 

5 

110  =  -2 

6 

111  =  -3 

7 

[ai  through  A(4-l-g^ 

is  requested. 

submitted  for  WORD  I. 

If  item  Cg(i,  j)  is  requested,  the  group  to  which  the  request  applies 
must  be  submitted  in  the  form  of  an  octal  word,  WORD  U.  WORD  II,  like  WORD  I. 


1. 


See  Edit  Input  sheet 
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is  first  represented  as  a  binary  number  and  then  converted  to  octal.  If  item 
Cg(i,  j)  is  not  requested.  +000000000000  must  be  submitted  for  WORD  II. 

Item  B(i,  j)  is  requested  by  specifying  the  plane  and  columns  for  which 
the  data  are  desired.  Input  data  consist  of  sets  of  three  consecutive  words 
(i,j’,j’’).  Each  set  selects  all  the  columns  in  plane  i  between  columns  j'  and 
j"  inclusive.  If  only  one  column  is  desired,  j"  =  j'  must  still  be  specified.  If 
two  or  more  sets  are  desired  for  a  single  plane,  the  sets  must  be  listed  in  order 
of  increasing  j  and  may  not  overlap.  Sets  must  be  listed  in  order  of  increasing 
plane  i.  The  last  set  must  be  followed  by  a  zero.  If  item  B(i,j)  is  not  request¬ 
ed,  a  zero  must  be  supplied  as  input. 

Item  Cg(i,  j)  is  requested  in  the  same  manner  as  item  B(i,  j).  The  sets 
of  (i,j',j")  follow  the  zero  after  the  B(i,  j)  sets.  The  last  set  must  be  followed 
by  a  zero.  If  item  Cg(i,  j)  is  not  requested,  a  zero  must  be  supplied  as  input. 
For  example,  if  neither  item  B(i,  j)  nor  item  Cg(i,  j)  is  requested,  the  third 
edit  request  input  card  contains  (.3252, 0, 0, ).  If  the  selected  columns  for  fluxes 
under  item  Cg(i,  j)  are  the  same  as  those  selected  under  item  B(i,  j),  a  one  may 
be  put  in  the  box  under  SEN  in  WORD  II,  and  a  zero  substituted  for  the  sets  of 
(i,  j',  j’)  's  for  item  Cg(i,  j). 


tiXiJni,  ti"  <)t  I’JDl  l'  PROGRAM  OUTPLj'I' 

TRIXY  ~ 

DATA  1  PLAMt  I  COLliMN  1  GROUP  RUN  1 

PROBLEM  NUMBER  1 CROjS  SECTION  DATA  TILE  NUMBER  ! 

K  POWER  DLliSlTY 

1  OOGOOOOO 

2  OOOOOui,'>) 

3  00000000 

POWER  INTEGRAL  00000000  POWER  DEf-lSIFY  00000000 


DATA  1  PLANc 
PROBLEM  NUMBER 


COLUMN  Z  GROUP  RUN  1 

10  CROSS  SECTION  DATA  FILE  NUMBER 


K  POWER  DENSITY 

1  28726732 

2  23390913 

3  11649406 

POWER  integral  3  20427648  POWER  DENSITY 


21789491 


DATA  1  PLANE 
PROBLEM  NUMBER 


1  COLUMN  3  GROUP  RUN  1 

10  CROSS  SECTION  DATA  FILE  NUMBER 


POWER  DENSITY 
1  16321599 
1  13479241 
1  10534715 


POWER  INTEGRAL 


4  12612843  POWER  DENSITY 


1  13453699 


DATA  1  PLANE 
PROBLEM  NUMBER 


2  COLUMN  1  GROUP  RUN  1 

R  10  CROSS  SECTION  DATA  FILE  NUMBER 
POWER  DENSITY 
28726732 
23390912 
11649405 


POWER  integral 


3  20427647  POWER  DENSITY 


21789490 


DATA  1  PLANE 
PROBLEM  NUMBER 


COLUMN  2  GROUP  RUN  1 

10  CROSS  SECTION  DATA  FILE 


POWER  DENSITY 
43367021 
35583599 
22139617 


POWER  INTEGRAL 


22423051  POWER  DENSITY 


NUMBER 


34168459 


DATA  1  PLANE 
PROBLEM  NUMBER 


POWtR  INTEGRAL 


3  COLUMN  3  GROUP  RUN  1 

R  10  CROSS  SECTION  DATA  FILE  NUMBER 
POWER  DENSITY 
1  16299149 
1  13569056 
1  10191957 

L  4  62846738  POWER  DENSITY  1  13407304 


I'AIA  4  PLANE 
PROBLEM  NllMHEK 


COLUMN  GI<OUP  I^UN  1 

10  lKOSS  SECTION  DATA  FILE  NUMBER 


NO  FISSION',  t'FF  SOURCE  NLUIRON 
{(OOOOOOO 
1  3  }  152  3'i 

32 1 3‘G  n 

4  5  4  /  4  4  (J  7 


Fiuurti  17  in 


EDIT  PROGRAM  OUTPUT  (con't.  ) 
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CATA  ::  PLAMf:  COLUMN  GROUP  1  RUN  1 

PROr^LLM  NUMPEK  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  FLUX  INTEGRAL 

1  6  23154808 

1001  6  11387804 

1002  4  45632656 

DATA  5  PLANE  COLUMN  GROUP  2  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  FLUX  integral 

1  6  17502494 

1001  5  93046648 

1002  4  32566869 

DATA  5  PLANE  COLUMN  GROUP  3  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  FLUX  INTEGRAL 

1  5  42182665 

1001  5  31505944 

1002  3  11430565 

DATA  6  PLANE  COLUMN  GROUP  1  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  AVERAGE  FLUX 

1  1  88208790 

1001  1  30987223 

1002  1  45632656 

DATA  6  PLANE  COLUMN  GROUP  2  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  AVERAGE  FLUX 

1  1  66676165 

1001  1  25318815 

1002  1  32566869 


DATA  6  PLANE  COLUMN  GROUP  3  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

INDEX  AVERAGE  FLUX 

1  1  16069587 

1001  85730459 

1002  11430565 

DATA  7  PLANE  1  COLUMN  2  GROUP  1  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

K  FLUX 

1  1  79071158 

2  1  63674468 

3  1  22435395 

4  00000000 


DATA  7  PLANE  1  COLUMN  3  GROUP  1  RUN  1 

PROBLEM  NUMBER  10  CROSS  SECTION  DATA  FILE  NUMBER  1 

K  FLUX 

1  2  15247520 

2  2  12542932 

3  1  50514745 

4  00000000 


Figure  17  (2) 
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mbol 

Mesh 

i 

j 

k 


i.  j,k 


I 

J 

K 

x(i) 

X(2) 

X(3) 


X.  GLOSSARY  OF  SYMBOLS 


Unit 


Possible 

Range 


Description 


integer  1  to  I 

integer  1  to  J 

integer  1  to  K 

integers 


integer 

integer 

integer 

cm 

cm 

cm 


Mesh  plane  perpendicular  to 

x(i) 

Mesh  plane  perpendicular  to 

X(2) 

Mesh  plane  perpendicular  to 


Location  of  mesh  point  at  inter - 
'  section  of  three  planes .  (Also 
used  to  identify  composition  and 
flux  boxes,  as  described  below). 
Maximum  i 
Maximum  j 
Maximum  k 
Cartesian  coordinates 
Cartesian  coordinates 
Cartesian  coordinates 


J 
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Symbol 

Mesh 

x<‘> 


x<'' 

J 


X?) 


M)  „(2)  (3) 

Xj  ,  Xj  ,  Xj^ 


^  (1) 


Unit 

cm 

cm 

cm 

cm 

cm 

cm 


A 


A 


(2) 

j 


(3) 

k 


cm 


cm 


Geometric 

\),x 


3 

cm 


Possible 

Range _ Description _ 

Coordinate  of  mesh  plane  i 
Coordinate  of  mesh  plane  j 
Coordinate  of  mesh  plane  k 
Coordinates  of  boundary  planes 

Coordinates  of  boundary  planes 

Spacing  between  mesh  planes 
i  and  i+1 

Spacing  between  mesh  planes 
j  and  j+1 

Spacing  between  mesh  planes 

k  and  k+1 .  Outside  the  mesh 

=  0,  where  /'  =  1, 
m 

2,  or  3,  andm'=0,  I,  J,  or  K. 

Volume  of  a  box  whose  corners 
are(i,j,k),  (i,j,k+l),  (i,j+l,k), 

(i,  j+1,  k+1),  (i+l,j,k),  (i+1,  j,  k+1). 


Ill 


Possible 

Symbol  Unit  Range 

Geometric 


V. 

1. 


j.lc 


3 

cm 


3 

cm 


Description 


(i+l,j+l,k),  (i+1,  j+l,k+l). 

This  box  must  be  of  uniform 
composition,  and  is  referred  to 
as  ’’composition  box”  (i,  j,  k). 

V.  u  ^  V.  is  the 

i,j,k  8  i,j,k 

octant  of  the  composition  box 

(i,  j,  k)  nearest  to  mesh  point 
(i,  j,  k).  This  elementary  volume, 
obtained  by  dividing  each  com¬ 
position  box  in  eight  equal  parts, 
is  referred  to  as  an  ’’octant” . 


V  =  V 

I.  j,  k  i,  0,  k 

V  =  V 

i.J.k  i,j,0 

V.  .  „  =  0). 

1,  J.K 


i.  j,  k 


I 


1 


I ,  m,  n=0 


i-  I,  j-m,  k-n 


is  the  total  volume  of  the  eight 


octants  adjacent  to  mesh  point 
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Symbol 

Geometric 


a 


(1) 

j.k 


a(2) 

i.k 


a 


(3) 
i.  j 


Possible 

Unit  Range 


Description 


2 

cm 


2 

cm 


2 

cm 


3 

cm 


(i,  j,  k).  This  volume  is  referred 


to  as  "flux  box”  (i,  j,  k). 


Area  of  the  projection  v. 

Ji  k 

perpendicular  to 

(1)  1  (2)  (3) 

“j.k'i 

(2) 

Projection  of  v.  X  , 

^ 

(2)  _  1  (1)  (3) 

“i,k~4  ^i 

(3) 

Projection  of  v  X 

J>  if  / 

(3)  1  (1)  (2) 

i,j“  4  ^i  ^j 

Each  of  these  areas  is  referred 


to  as  a  "quadrant",  since  four 


quadrants  comprise  one  face  of 


a  flux  box. 

(  I ') 

a  .  ,  =  0,  where  m'  or 

m^,n 

n'  =  0.  I,  J.  K, 


Total  volume  of  regions  con 


taining  fissionable  material 
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Symbol 

Geometric 

•  I 

V 


F 

Z 


Indexes 

g 

f 

g 


Possible 

Unit  Range 


Description 


integer  1  to  G 

integer  1  to  G 


I.  J.  K  .. 

V 


Fissionable  volume  in  flux  box 
(i.  j.  k), 


^  V. 

,  m,  n=0,  1  ^ 


,  j-m,  k-n 


i-  , j-m, k-n’ 

Summation  to  be  carried  out  only 

over  that  portion  of  the  mesh 

containing  fissionable  material. 

F  F 

p  i  max  j  max 

n  =  I,  Zp 

i=i  min  j=j  min 
,F 

k  max 
k=k‘ min 


R®  refers  R  to  energy  group  g. 
In  all  cases,  fission  cross 

t 

section  is  zero  for  g  <  g  . 
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Symbol 

Unit 

Possible 

Range 

Description 

Indexes 

G 

integer 

1  to  1 5 

Maximum  g. 

(t) 

integer 

0  to  63 

^  refers  R  to  the  t^^  mesh 

(T) 

integer 

0  to  63 

sweep  during  a  group  iteration. 

Maximum  t  (the  last  mesh 

(r) 

integer 

0  to  63 

sweep  of  a  group  iteration). 

( r ) 

R'  refers  R  to  run  r. 

^ ii j.  k 

integer 

0  to  63 

Composition  index  (see  symbols 

q 

integer 

0  to  7 

under  "Composition"). 

Identification  of  an  octant  in  a 

flux  box.  The  nximber  q  is  the 
decimal  equivalent  of  the  binary 


number  (1,  m,  n). 
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Symbol 

Unit 

Possible 

Range 

Description 

Indexes 

1 

integer 

0  to  1 

Thei,  m,  or  n=0  denotes  the 

m 

integer 

0  to  1 

positive  direction  along 

or  X^^^  respectively,  and 

m,  or  n=l  denotes  the  negative 

n 

integer 

0  to  1 

direction.  Positive  and  negative, 

in  this  case,  are  referred  to  an 

origin  at  mesh  point  (i,  j,  k) , 

Composition 

1. 

-1 

cm 

Moderating  (slowing-down) 

-1 

cm 

cross  section. 

Fission  cross  section. 

1, 

cm 

Diffusion  coefficient. 

1,  J,K 

A®.  . 

-1 

cm 

Removal  cross  section  =  M  + 

1.  J,K 

1 

i,  j,k 

integer 

0  to  63 

absorption  cross  section. 

Composition  index  of  composition 

box  (i,  j,  k).  For  each  index,  G 

sets  of  (M®,  F®,  D^,  and  A®) 
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Symbol 

Composition 


Possible 
Unit  Range 


i.  j.k 


integer  0  or  1 


7.(l)g 

7(l)g 

^(2)g 

^(2)g 

v(3)g 

^(3)g 


Description 


must  be  specified.  If  composition 
box  (i,  j,  k)  falls  outside  of  the 
mesh  (i,  j,  k  <  1  or  >  I-l,  J-1, 
K-1),  then  each  M,  F,  D,  A  io 
zero. 

If  composition  box  (i,  j,  k)  contains 
fissionable  material  use  1;  other¬ 
wise  use  0 . 

Ratio  of  outward  current  to  flux 
at  a  boundary  plane  (applicable 
when  a  current  boundary  condi- 


S  D  Of 

tion  is  specified),  Z®( - 

a  property  of  the  medium  out¬ 
side  the  boundary  plane,  which 
is  assumed  constant  over  the 


is 


entire  plane . 
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Symbol 

Coefficients 

p(l)g 

i.j.k 


Unit 


cm 


Description 


c“>S  =  d8  ,  . 

i.j.k  i,j-m,k-n  j-m,k-n/  i 

m, n=u  ' 


,(2)g 

l.j.k 

,(3)g 

Xj.k 


cm 


cm 


c(2)g  D«  ^(2)  /  (2) 


c(2)«  -  ^  dS 

i-^.j-m.k  i-tj 


(3)  /  i3) 

m  /  k 


E 


g 

i.j.k 


cm 


'S  ^  S  V 

fCrT  n  i-i.j-m,k-n  i-i.j-m,k-n 
jf,m,n=U 


I.j.k 


cm^  B®  .  ,  = 

i.j.k  I.j.k  i'-l,j,k  I.j.k  i,j-l,k 

+  u  ,  +  E.  .  , 

I.j.k  i.j.k-1  i.j.k 


1,  J.k 


1 


g 


cm  .  ,  =  'y  .  1,  '  •  0  ^  rv.  ir  T. 

I.j.k  i-jl,j-m,k-n 


i.j.k 


cm 


1 


p  g  =  ^  pg  V 

i.  ji-k  „  - n  i--*.j-m.k-n  i-J?,  j-m,  k-n 

X,  Tci,  n=u 


Boundary  Coefficients 


.(l)g 

"i'.  j.k 


cm 


(l)g  ^  ^(Ijg 
i’.j.k  i' 


m,  n=0 


j-m,  k-n  i'=0  or  I 


118 


Symbol 

Unit 

Boundary  Coefficients 

c(2)g 

i.j'.k 

2 

cm 

c(2)g 

i.j’.k 

c(3)g 

i.j.k' 

2 

cm 

c(3)g 

i.j.k' 

Description 


a 


H,  m=0 


(2) 

'i-i,k-n  j' 

^(3) 

i-i,j-m  k 


Components  of 
Coefficients 

d(i>g 

i.  j.k 


P)g 

i. 

l(3)g 

i.  j,k 


rr,^ 

j,k 


(l)g 
i,  j,k 


J2)g 
i,  j.k 


cm 


cm 


cm 


cm 


cm 


cm 


cm 


=  D®  A^^^  /  4  A^^^ 

i.j.k  i,j,k  j  k  /  i 


d(2)g  _  Qg 

i.j.k  i.j.k 


4’/ 4 

1  k 


<3)«  =D«  A<1)  A  (2) 


a'.^> 

J 


(3) 


a 


i,j,k  i,j,k  i 

«  =A« 


A';'  A'“'/4  A 
j  /  k 


i,j,k  i,j,k 
m®  .  .  =  M® 


a‘'>  A<^)/8 
1  J  k  / 


i.j.k  i.j.k  “i  ''j 

i.j.k  i.j.k  i.j-l,k 

r<2>g  =d<2)g  +d<2)g 

i.j.k  i.j.k  i-l,j,k 


A<»  A<^> 


i.j.k 


cm 


5«  =  5“  aS 

ifs^o 


=  0  or  J 


'  =  0orK 
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Symbol  Unit 


Components  of 
Coefficients 

i.J.k 


2 

cm 


Flux 


i.  j.k 


2 

cm 


neutrons 

2 

sec -cm 


Source  and  Loss 


i, 


neutrons 

sec 


Description 


i.j.k 


1,  J.k 


1 

"IZ 

i,m=0 

i.  j.k 


i-i,  j-m,k 

i,j,k  i,j,k  i,j,k 


Flux  in  group  g  at  mesh  point  i,  j,  k  after  it 
has  been  improved  t  times  during  the  group 
iteration.  is  the  flux  estimate  at 

i,J,k 

the  beginning  of  the  group  iteration . 

-  J  - 

‘Po.j.k  ‘^H-l,j,k“  '^i,  0,k'  ^i,J+l,k 

Source  of  fission  neutrons  in  flux  box  (i,  j,  k). 

^  ,  is  normalized  to  an  average  of  one 

i,j,k 

neutron/sec  over  all  fissionable  regions . 

7 

Lj»k  (i,  j,k)q 
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Symbol  Unit 

Source  and  Loss 


S* 

i,  j.k 


neutrons 

sec 


%S* 
i,  j.k 


neutrons 

sec 


a** 

^  •  1, 

l.J.k 


neutrons 

sec 


i,j,k 


neutrons 

sec 


t) 

i,  j.k 


neutrons 

sec 


i,  j,k 


neutrons 

sec 


og(r) 

i,  j,k 


neutrons 

sec 


Description 


Fraction  of  the  fission  neutrons  that  are  born 
in  group  g.  0  <  <  1- 

Fission  source  at  point  (i,  j,  k),  not 
normalized. 

* 


Contribution  of  group  g  to  ife  .  .  .  . 

ii  Ji  -If 

,S*  _  tj'g 


i.j.k  i,j,k  i,j,k 


i.j.k 


g=g 


S®* 

i.J.k 


Neutron  source  in  flux  box  (i,  j,  k),  in  group 
g,  due  to  both  fission  and  degradation . 

S®  .  .  =  ¥®  .  .  + 

i,j,k  Jt.  i,j,k  1,  j,k  1,  J,  k 

*  =  S  /  f 

i.j.k  i,j,k/ 

Estimates  of  the  above  quantities  after 
completion  of  run 

Estimates  of  the  above  quantities  after 


completion  of  run 
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Symbol 

Source  and  Loss 

■  t 

i,j,k 


Unit 


Description 


u(  S®)  neutrons 


(l)g  (l)g 
^1  ’ 


T  (2)g  .  (2)g 

^1  ’  h 


,  (3)g  .  (3)g 

^1  ’ 


sec 


Neutron  removal  and  leakage  from  flux  box 
(ij.k),  in  group  g. 


=  bS  .  ,  , 

ifj.k  i,j,k 

1,  Jflc 

+  c<^>« 

i+l,j,k  i-l,j,k 

‘^i-1,  j.k 

l.JflC 

,g  .  p(2)g 

“^ij+l.k  ^  '"ifj-lfk 

,  V 

^ifj-lflc 

i.J.lt 

,g  ,  p(3)g 

‘^i.j.k+l  i,j,k-l 

refers  L®  .  ,  to  mesh  sweep 

ifjfk  ^ 

t  during  a  group  iterations . 

^  ®  S®)  ^  L®  .  .  (  0)  -  S® 

if  j»  k  i,  j,  k  i,  j,  k 

Neutron  leakage  throu^  boundary  plane 

For  specified  current  condition: 

I  (i)  _  p(i)g 

^  j^l 


=  ,  ate. 

‘ 
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Syml?ol  Unit 

Source  and  Loss 


neutrons 

sec 


E)escription 


Neutron  leakage  through  boundary  plane 
For  specified  flux  condition: 


(i)g_4^  ^(i)g 


U-'==  >  C',"?  ,  ,.(  <-?  ,  .  .  J,  etc, 


-I  H.J.k  I, 

Neutron  leakage  from  the  mesh  in  group  g. 


g 


Power 

g 


fissions  Power  density  in  quadrant  q  of  flux  box 

3 

SCC^CIU 

(i,  j,  k)  due  to  fissions  induced  by  group  g 


neutrons . 


fissions 


sec- cm 


P®  ,  =  P*  o* 

‘^<i.j.k)q  i-l,j-ni,k-n  l.j.k 

g 

Power  density  ^jf(i.j,k). 


neutrons 


V 


fission 


Symbol  Unit 

Miscellaneous  Factors 


S 

P 


? 


£ 


Rg(t-1) 


Description 


Flux  change  criterion  for  terminating  flux 


iteration 


g(t-lh2 

V 


< 


g 


Power  change  criterion  for  final  run 


I,J.K 
i,  j,k=l 


(r>  _  (r-1)  2 

'^i,j,k)q  (i,j,k)q' 


<  O 


Criticality  error  criterion  for  final  run. 

Renormalization  factor  after  (t-1)^  mesh 


sweep . 


a 


Liebmann  Acceleration  Factor 
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Symbol  Unit 

Miscellaneous  Factors 


e 


(r) 

P 


e 


(r) 


c 


1 


t 

max 


Description 


Flux  change. 


Power  change, 

'P  OTTSl  ^ 


Source  change  between  successive  runs, 
F 


(r) 


(r) 

'’i.j.k 


(r-1) 

j.k 


Convergence  criterion  for  terminating  the 


runs, 

(r>  ^  . 

e  ^  • 

Percentage  of  initial  total  improvement  to 
serve  as  a  secondary  convergence  criterion 
for  terminating  the  flux  iterations . 

Maximum  number  of  sweeps  allowable  in  the 
flux  iterations  (tertiary  convergence  crit¬ 


erion)  . 
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Symbol  Unit 

Miscellaneous  Factors 


Programming  Conventions 
r  integer 

J'  integer 

integer 

integer 

integer 


K' 


g 


Description 


Frequency  print  regulator  and  frequency 


regulator  for  computing  error  ratio. 


(t> 


(t-Cj) 


in  flux  iterations . 

Percentage  of  the  fissionable  volume  by 
which  the  computed  volume  and  the  estimated 
volume  may  differ .  If  the  difference  is 
greater  than  this  the  machine  stops . 


Number  of  planes  on  which  flux  is  computed. 
Number  of  columns  per  plane  on  which  flux 
is  computed. 

Number  of  points  per  column  on  which  flux 
is  computed. 

Group  number  of  the  flux  guesses  to  be 

improved  by  the  Group  Sweep  Program  . 

P  i 
g  =  g  +  1. 
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Symbol 

Programming 

C 


RCW 


Unit 

Conventions 

integer 


binary 

word 


Description 


0  if  K  <36 
1  if  K  >36 

(C  indicates  whether  the  sentinels  for  each 
column  are  carried  as  one  or  two  words . ) 
Reactor  Control  Word:  a  12  bit  representation 
of  boundary  conditions  on  the  six  faces  of  the 
reactor . 

odd  bits  denote  0  -  flux  boundary 

1  -  current  boundary 
even  bits  denote  0  -  zero  boundary 

1  -  non -zero  boundary 


35, 

34  refer  to 

i 

=  1 

33. 

32 

j 

=  1 

31, 

30 

k 

=  1 

29, 

28 

i 

=  I 

27, 

26 

j 

=  J 

25. 

24 

k 

=  K 
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Symbol  Unit  _ Description _ 

Programming  Conventions 

SFGW  binary  Source  Flux  Guess  Word:  It  is  zero  if  a 

word 

rough  power  guess  is  supplied,  or  non-zero 
if  the  results  of  a  previous  run  are  to  be  used. 
CC  integer  Convergence  Code  for  Group  Sweep  Program. 

1  -  primary  convergence  test  met 
CC  =  2  -  secondary  convergence  test  met 
3  -  tertiary  convergence  test  met 
A  block  of  14  words  of  control  input  are  stored  in  memory  and  preserved 


there  by  all  succeeding  programs .  This  block  contains: 


Decimal 

Location 

Decrement 

Tag 

Address 

29 

a  (floating  point  form) 

30 

g’ 

0 

G 

31 

gP' 

4 

RCW 

32 

.F 

1  max 

0 

r 

33 

.F 

1  mm 

0 

J' 

34 

.F 

j  max 

0 

K' 

g^  =  1  at  the  beginning  of  each  run.  g^  is  increased  by  1  at  the  end 
of  a  group  iteration  by  the  Group  Sweep  Program  before  control  is  transferred 
to  the  Group  Source  Program . 


1Z8 


Decimal 

Location 

35 

36 

37 

38 

39 

40 

41 

42 


Decrement  Tag 

F 

j  min  0 

F 

k  max  0 

F 

k  mm  0 

1 

Tape  A  0 

Tape  C  0 

Tape  E  0 

0  4 


Address 

I 

J 

K 

Tape  B 
Tape  D 

2 

Flux  Guess  Tape 
Tape  F 

1  (Hrji.4) 


^  Master  tape 
2 

The  Flux  Guess  Tape  is  either  tape  B  (and  the  sign  of  40  is  -)  or  tape  E 
(and  the  sign  of  40  is  +) . 
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APPENDIX  A:  INPUT 

INSTRUCTIONS  FOR  FILLING  OUT  INPUT  FORMS 
General  Instructions 

The  lines  on  most  of  the  input  sheets  are  divided  into  two  parts .  The 
upper  half  of  a  line  contains  identification  and  headings  which  serve  as  guides  to 
the  person  preparing  the  input .  The  lower  half  contains  the  numerical  data 
which  are  punched  by  the  keypunch  operator .  K  no  letter  or  number  is  given 
in  column  one,  that  column  is  to  be  left  blank  unless  this  is  the  last  card  of  a 
record.  An  ei^t  in  column  one  indicates  this  is  the  last  card  of  a  record.  This 
card  may  also  have  other  punches  in  column  one  and  must  contain  data.  There 
can  be  no  blanks  within  a  record,  if  a  value  is  zero,  enter  the  number  0.  Except 
for  RPG  (Rough  Power  Guess),  no  records  can  be  omitted.  The  maximun  number 
of  digits  is  eight  plus  exponent .  Commas  cannot  be  used  to  separate  the  digits 
in  a  number .  Floating  point  numbers  can  be  entered  in  several  ways .  For 
example,  jr  can  be  entered  as: 

3.1415927 

+31416/-4 

.03142/2 

Fixed  Point  Geometrical  and  Physical  Data 

Memory  Size  -  Must  be  given  to  the  nearest  thousands,  i.e.,  8000, 


16,000,  32,000. 
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CSFN  -  Cross  Section  File  Number,  The  file  of  cross  section  records 
to  be  selected  from  the  configuration  and  cross  section  tape , 
RCW  -  Reactor  Control  Word.  This  word  is  computed  by  defining  the 
boundary  conditions  vrtiich  are: 

0  -  zero  flux 

1  -  zero  current 

2  -  non  zero  flux  (not  coded  at  present) 

3  -  non  zero  current 

The  appropriate  values  are  multiplied  by : 

1  for  i  =  1 

4  for  j  =  1 
16  for  k  =  1 
64  for  i  =  I 

256  for  j  =  J 
1024  for  k  =  K 

The  sum  of  these  products  is  the  RCW . 

SFGW  -  Source  Flux  Guess  Word.  Positive  for  elastic  scattering  and 
negative  for  inelastic  scattering  (not  coded  at  present) ,  The 
value  entered  is  zero  if  a  flux  guess  is  to  be  supplied  from  a 
previous  problem.  It  is  one  if  a  rough  power  guess  is  supplied 
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Floating  Point  Input 

V„  -  Fissionable  volume.  If  an  estimate  is  not  made  this  should  be  0. 

F 

%V  -  This  is  the  permissible  error  in  V  and  is  a  decimal  fraction. 

F  F 

List  A  and  List  B 

Both  of  these  lists  must  always  be  supplied.  If  one  of  the  lists  is  to  be 

disregarded,  the  number  100  must  be  written  as  the  first  word  of  the  record. 

Then,  column  1  has  an  8  punch  in  addition  to  the  letter  D.  The  column  headed 

with  an  asterisk  is  to  be  used  only  if  100  is  written  in  that  space. 

Coarse  Mesh  and  Rough  Power  Guess 

The  coarse  mesh  supplied  for  the  rough  power  guess  must  not  be  larger 

than  the  actual  mesh ,  The  first  coordinates  of  the  rough  power  guess  cannot  be 
F 

greater  than  K  there  are  less  than  six  values  given  for  any  of  the  coordin 

ates,  the  remaining  blanks  on  that  line  should  be  disregarded  by  both  the  person 
preparing  the  input  and  the  keypunch  operator. 

Configurations 

These  must  be  listed  in  order  of  increasing  C  starting  with  C  =  1 . 

n  n 

ORDER  OF  INPUT 

First  Record^;  Fixed  Point  Geometrical  and  Physical  Data 

Current  Boundary  Conditions 
Floating  Point  Input 

^  If  a  value  other  than  2.46  is  desired  for  V  ,  it  can  be  entered  with  the 
first  record,  e.g. ,  if  =  1 .5  the  card  would  be  punched  Cl,  1 . 5, 
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Second  Record;  List  A 


Third  Record;  List  B 


Fourth  Record;  Mesh  Spacings  and  Coordinates 
Fifth  Record  Edit  hiput 

Sixth  Record;  Rou^h  Power  Guess  (if  supplied) 

OMEGA  (one  card  for  each  run  -  in  card  reader) 

Two  blank  cards 

INSTRUCTIONS  TO  MACHINE  OPERATORS 


It  is  necessary  that  the  operator  know; 

1 .  If  input  is  on  cards  or  on  tape, 

2.  If  a  control  restart  deck  is  to  be  punched, 

3.  If  the  Master  Tape  is  to  be  duplicated,  and 

4.  If  Monitor  printing  is  wanted  during  the  Group  Sweep.  At  present.  Omega 
must  be  read  in  from  cards .  Also,  if  Vp  is  supplied,  a  tolerance  must 
be  given .  If  the  Vp  computed  is  outside  this  tolerance,  the  machine 
stops  and  this  difference  is  printed.  If  this  difference  is  small  and  it  is 
decided  to  continue,  this  must  also  be  indicated. 


TRIXY  INPUT  DATA 
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(100,)  IMMEDIATELY  FOLLOWS  LAST  j,  REGARDLESS  OF  COLUMN 


TRIXY  INPUT  DATA 
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AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  THIS  INPUT 


u 

(0 

2 

1 — I — 1  1 — 1 — 1  1 — i~r'i 

139 

1 — m  1 — 1 — 1  1  -r.i 

M 

vM 

HHH 

4H 

^■Mp 

Mfl 

^■Hp 

^Wp 

Hi 

mm 

^Bp 

<w 

hhum^^H 

TiUXY  INPUT  DATA 
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Column  1  of  the  first  two  cards  contains 


TRKY  INPUT  DATA 
Configurations 


142 


Column  1  will  be  left  blank  on  all  continued  lines  EXCEPT  the  last  line  of  a  config^ation,  then  it  will  have  an 
written  in  it. 

The  first  word  of  the  configuration  can  be  either  a  fill  (no  sign)  or  a  skip  (-  sign). 

100  must  be  written  at  the  end  of  each  column  and  1000  at  the  end  of  each  configuration. 


TRIXY  INPUT  DATA 
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AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  EACH  RECORD. 


TRIXY  INPUT  DATA 
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AN  8  MUST  BE  IN  COLUMN  1  of  LAST  DATA  CARD  OF  EACH  RECORD. 


Page  1  of _ NF 
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U 

2 

2 

w 

53 

tX4 

I 

§ 

2 


Q 


§ 

cd 

H 


C*5 

T3 

n 

o 

o 

(D 

Oi 


«>> 

C 

o 


o 

<u 

CO 

CO 

CO 

o 

u 


AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  EACH  RECORD 
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Z 


AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  EACH  RECORD. 


Page  1  of _ CBC 
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AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  EACH  RECORD 


Page _ of  _ CBC 
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AN  8  MUST  BE  IN  COLUMN  1  OF  LAST  DATA  CARD  OF  EACH  RECORD. 


APPENDIX  B- 


DECK  MAKE-UP 

GE  BILD  1 

GE  BILD  2 

PKl 

(75  cards) 

LTW 

(12  cards) 

TRXSLOAD 

(1  card) 

PKl 

(75  cards) 

MT-1 

(71  cards) 

Blank  card 

G.SOR  (86  cards) 


Blank  card 

SWEEP  (92  cards) 
Blank  card 

MT-3  (67  cards) 

FISS  (73  cards) 

Blank  card 

EDIT  (139  cards) 

Blank  card 

MT-2  (110  cards) 

4  blank  cards 


GENERATION  OF  LIBRARY  TAPE 


(3000  words) 

(2878  words) 

(1825  words) 

(1559  words) 
(1605  words) 

(3079  words) 

(2728  words) 
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OPERATING  INSTRUCTIONS 

Place  blank  tape  on  unit  3,  Sense  switches  2  and  3  must  be  down.  CLEAR 
and  LOAD  CARDS;  final  stop  (33l6)  . 

O 

DESCRIPTION 

Errors  in  reading  cards  or  writing  tape  are  indicated  on  the  printer.  As 
each  program  is  written  on  tape,  the  number  of  words  in  the  program  is  printed 
out.  If  this  does  not  check  with  the  word  count  given  above,  an  error  in  deck 
make-up  is  indicated.  This  program  loads  both  SAP  and  CAGE  binary  cards  It 
does  not  continue  if  a  card  with  an  incorrect  check  sum  is  discovered. 


151 


APPENDIX  C:  GENERATION  OF  CONFIGURATION  AND  CROSS  SECTION 

DATA  TAPE 

DECK  MAKE-UP 

Binary  Cards 

1 .  BILD  I. 

2.  BILD  n. 

3.  GEPRO  (71  cards).  This  is  a  special  copy  of  GEPRO. 

4.  Tape  4  (21  cards). 

5.  Transfer  card  (3173)^. 

Decimal  Cards 

1 .  Configuration  data  cards.  One  configuration  per  record. 

2.  End  of  File  card  (a  blank  card  with  an  8  in  column  1). 

3.  Cross  Section  Data  cards  for  File  I. 

4.  End  of  File  card. 

5.  Subsequent  Cross  Section  Data  cards  for  respective  files .  Each 
file  followed  by  an  EOF  card. 

6.  Two  End  of  File  cards  follow  the  last  Cross  Section  file. 
OPERATING  INSTRUCTIONS 


If  the  data  cards  are  put  on  BCD  tape,  the  tape  must  be  put  on  unit  2. 
Regardless  of  whether  tapes  or  cards  are  used,  the  output  will  be  on  tape  unit  4 . 
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Sense  switch  2  must  be  UP  if  data  are  on  tape,  DOWN  if  on  cards ,  Sense  switch 
3  must  be  EXDWN. 

The  following  errors  are  detected  by  this  Deck: 

1.  Incorrect  number  of  words  in  a  cross  section  record. 

2.  Incorrect  order  of  cross  section  files. 

3.  Cross  section  A  equals  zero. 

4.  Incorrect  placement  or  omission  of  End  of  File  card, 

5.  The  number  of  words  in  record  2  is  not  a  multiple  of  4G. 

6.  The  number  of  words  in  record  3  is  not  a  multiple  of  3G. 

7.  The  number  of  words  in  record  4  is  not  a  multiple  of  G. 

8.  Failure  of  tape  4  to  write  correctly. 

9 .  Any  mispunched  data  cards . 

If  any  of  these  errors  occur,  corrections  must  be  made  and  then  the  program 

re-loaded.  Final  stop  (3446)  , 

8 
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APPENDIX  D:  PREPARATION  OF  A  NEW  FLUX  GUESS  TAPE 


DESCRIPTION 


This  program  to  prepare  a  new  flux  guess  tape  has  as  input  tape  E  from 
a  previous  run  and  record  1  of  the  input  for  a  new  case .  It  generates  a  new  flux 
guess  tape  with  the  history  records  corrected  to  correspond  to  the  new  input. 
The  f  file  is  omitted  from  the  new  tape.  Only  G  files  of  flux  are  copied  onto 
the  new  tape .  Since  the  files  containing  the  rouj^  source  guess  and  flux  are 
copied  directly,  the  following  words  of  input  for  the  two  cases  must  be  the 


same: 

*•  J’  ‘L'  Cax’  >L'  iLc'  “Ln-  ''max’ 

The  G  for  the  new  case  cannot  be  greater  than  the  G  for  a  previous  case. 

DECK  MAKE-UP 

Binary  cards: 


BILD  I 
BILD  II 

PKl  (75  cards). 

Tape  5  (12  cards). 

Decimal  cards: 

Record  one  of  the  input  for  the  new  case  (16  cards) 
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OPERATING  INSTRUCTIONS 

The  results  from  the  previous  case  must  be  placed  on  tape  6 .  The  new 

tape  is  generated  on  tape  5.  Sense  switches  2  and  3  must  be  down.  Any  stop 

except  the  final  stop  (3334)  ,  indicates  either  machine  error  or  errors  in  the 

8 

decimal  cards.  There  should  be  a  print-out  which  indicates  the  type  of  error. 
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APPENDIX  E  OPE R  >TING  MANUAL 


INTRODUCTION 

The  purpose  of  this  manual  is  to  provide  operating  instructions  for  the 
TRIXY  program.  It  is  suggested  that  this  operating  manual  be  kqjt  in  a  loose 
leaf  notebook  for  convenient  reference  while  operating  the  program , 

The  operation  of  TRIXY  differs  from  standard  operating  procedures  used 
by  most  installations .  Hence,  it  is  recommended  that  all  persons  concerned 
with  programming,  preparation  of  input,  and  machine  operation  of  TRIXY  read 
this  description. 

GENERAL  DESCRIPTION 

TRIXY  is  designed  to  run  for  several  hours  Therefore,  it  is  important 
that  the  operating  instructions  and  the  instructions  from  the  customers  be  followed 
carefully.  All  tape  units  must  be  cleaned  in  preparation  for  running,  and  only 
full  reels  of  tape  should  be  used.  Comments  are  printed  to  advise  the  operator 
of  the  progress  of  the  program  and  to  provide  instructions  whenever  intervention 
by  the  operator  is  required. 

The  TRIXY  Library  Tape  should  be  placed  on  unit  1  and  the  configuration 


and  Cross  Section  Tape  on  unit  4 ,  If  the  input  is  on  (BCD)  tape,  this  should  be 
placed  on  unit  5.  Blank  tapes  should  be  placed  on  units  2,  3  and  6,  If  the  input 
is  on  cards,  die  deck  should  be  placed  in  the  card  reader  and  sense  switch  2 


depressed ,  "OMEGA”  data  should  be  placed  in  the  card  reader  behind  the  input , 

If  input  is  on  tape,  the  OMEGA  cards  should  be  placed  directly  in  the  card  reader . 
If  the  customer  desires  a  Control  Restart  deck,  depress  sense  switch  4,  The 
program  stops  after  the  deck  is  punched  to  allow  the  control  rod  positioning  to  be 
checked.  This  deck  can  also  be  used  to  restart  the  master  tape  generation.  If 
the  deck  is  to  be  used  for  restart  only,  then  the  program  should  be  continued  by 
pushing  START. 

After  the  input  has  been  read,  the  program  reads  from  tape  4  and  writes 
on  tape  3.  Several  passes  are  made  to  sort  from  tape  3  to  tape  2.  This  operation 
may  take  more  than  an  hour .  When  this  section  is  finished,  the  control  restart 
deck  is  punched  if  requested  Information  is  transferred  from  tape  1  to  5.  The 
machine  stops  and  tapes  1  and  4  must  be  replaced  as  indicated  by  instructions  on 
the  printer.  At  this  point,  specific  sense  switch  settings  must  be  made  to  control 
the  operation  of  the  program  ,  It  is  recommended  the  sense  switches  be  set  so 
that  a  duplicate  of  the  master  tape  is  generated  and  monitor  printing  is  provided 
during  the  Group  Sweep .  (Group  Sweep  printing  does  not  take  extra  time,  except 
for  extremely  small  cases  (IJK  <  1000)  Printing  is  done  while  a  tape  is  rewound.) 

The  generation  of  the  master  tape  is  performed  Data  is  computed  and 
written  on  tapes  1,  3,  4  and  6  (if  a  duplicate  of  the  master  tape  is  being  written ,) 
The  group  number  is  printed  and  data  transferred  from  tapes  1  and  5  to  3  (and 
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possibly  6).  This  process  is  repeated  for  each  group  At  the  end  of  the  master 
tape  generation,  data  are  transferred  from  tapes  1,  2,  4  and  5  to  3  (possibly  6), 

If  a  rough  source  guess  is  required,  it  is  computed,  and  tapes  4  and  5  are 
used.  This  computation  is  performed  while  rewinding  the  master  tape  (s)  The 
machine  stops  If  only  one  master  tape  has  been  generated,  it  is  possible  to 
provide  a  copy  of  the  master  tape  by  pushing  START.  If  a  copy  is  not  wanted, 
put  sense  switch  1  up  and  push  START. 

The  Group  Source  Program  is  read  from  the  master  tape  and  control  is 
transferred  to  this  program,.  In  case  a  restart  procedure  is  being  executed,  the 
Group  Source  Program  positions  the  tapes  properly,  reads  in  the  Group  Sweep 
Program,  and  transfers  control.  If  a  restart  procedure  is  not  being  executed, 
values  of  the  scattering  coefficient  (H)  for  one  plane  are  read  into  memory  from 
the  master  tape  (tape  A),  A  plane  of  iterated  flux  (  f/>)  is  read  from  tape  B,  and 
the  values  of  'V  are  transferred  to  tape  E  for  the  Group  Sweep  Program.  H0 
is  calculated  and  stored  on  tape  C.  This  process  is  repeated  for  all  planes 
(This  first  step  is  omitted  for  the  first  group.) 

A  plane  of  fission  source  (  S  )  is  read  from  tape  E .  The  values  are 

O' 

multiplied  by  X  and.  for  groups  other  than  the  first,  are  added  to  the 
corresponding  value  of  H<P  from  tape  C  to  determine  the  group  source  (S). 

Blocks  of  ErrC  coefficients  are  read  in  from  tape  A,  combined  with  values  of 
S,  and  written  on  tape  D  as  blocks  of  SErrC’s .  If  a  flux  guess  is  not  available 
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for  the  succeeding  group,  ^  is  calculated  from  S  and  E  and  written  on  tapes 
B  and  E .  This  process  is  repeated  for  all  planes . 

The  Group  Sweep  Program  is  read  from  the  master  tape  and  control  is 
transferred  to  it . 

At  the  beginning  of  the  Group  Sweep  Program,  one  line  is  printed.  This 
line  is  printed  at  the  end  of  each  sweep  and  again  MAien  convergence  is  reached. 
During  each  sweep,  values  for  the  flux  ( <i)  axe  read  from  either  tape  B,  C, 
or  P  and  written  on  one  of  these  tapes .  Values  of  SErrC  are  read  from  tape  D. 
At  the  end  of  each  sweep,  tapes  B,  C,  D,  and  F  are  rewound  while  the  line  is 
printed.  When  convergence  is  reached,  the  Group  Source  Program  is  read  from 
the  master  tape. 

Control  alternates  between  Group  Source  and  Group  Sweep  until  the 
Fission  Source  Program  is  reached.  Fission  Source  reads  data  from  the  master 
tape,  tape  E  (the  tape  containing  the  flux  and  fission  source  guesses),  and  a  flux 
tape  prepared  by  Group  Sweep .  One  "OMEGA”  card  is  read  in  from  the  card 
reader .  A  new  tape  E  is  prepared .  Information  concerning  the  convergence 
and  the  tape  number  of  the  new  tape  E  is  printed,  and  the  machine  stops .  At 
this  point,  control  is  turned  over  to  Group  Source,  unless  the  fission  source  has 
converged  or  an  Edit  is  requested  by  depressing  switch  6. 

If  control  is  turned  over  to  the  Edit  Program  the  tape  to  be  used  for 
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output  is  indicated  on  the  printer .  This  program  first  sorts  the  computed 
results  from  tape  E  and  the  master  tape  onto  the  working  tapes;  then  it  writes 
the  results  on  the  output  tape.  This  output  may  take  more  than  one  reel.  K 
more  than  one  reel  is  needed,  the  machine  should  be  put  on  manual  before  the 
end  of  the  tape  is  reached.  This  causes  the  machine  to  stop  (not  a  programmed 
stop) .  An  end  of  file  should  be  written  on  the  output  tape,  another  output  tape 
turned  on,  and  the  machine  started  again.  The  end  of  TRIXY  is  indicated  on  the 
printer . 

NORMAL  OPERATING  PROCEDURE 
1  Arrange  the  tapes  as  follows: 

Unit  Description 

1^  TRIXY  Library  Tape 

2  blank 

3  blank  (labeled  master  tape) 

4^  Configuration  and  Cross  Section  Tape 

5  (BCD)  input  or  blank 

6  blank  (duplicate  of  master  tape) 

2.  Set  sense  switches  as  follows: 

1  up 

^This  tape  must  have  the  rite  protect  ring  removed. 
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(Set  sense  switches  continued) 

2  up  input  on  tape  5 

down  input  in  card  reader 

3  up 

4  up  no  restart  deck 

down  a  Control  restart  deck  is  requested 

5  up 

6  up 

3 .  Place  the  input  in  the  card  reader  (if  it  is  not  on  tape) . 

4  „  Place  the  OMEGA  deck  in  the  card  reader . 

5  Push  CLEAR  and  LOAD  TAPE 

Following  is  a  listing  of  the  printouts  and  stops  which  occur  during  a  normal  run 
of  TRDCY . 

Master  Tape  Generation 

Stop  Printout 

None  "The  I  Planes  of  'Composition  Indices  are  Finished’ . " 

161  "The  Tapes  have  been  set  up  —  " 

At  this  point  switches  2  and  3  only  should  be  depressed 
unless  otherwise  instructed 

None  The  group  number  is  printed  out  as  each  group  is  finished 


161 


Stop  Printout 

I6l  "A  copy  of  the  master  tape  is  on  unit  6.  Replace  and  save 

tape  6  " 

It  is  also  recommended  that  the  file  protect  ring  be  removed 
from  tapes  3  and  6. 

None  "Remove  Deck  from  Punch.  - " 

Group  Source  Program 

No  stops  or  printouts 
Group  Sweep  Program 

A  line  of  information  is  printed  at  the  beginning  of  the  Group  Sweep 
Program.  At  the  end  of  each  sweep,  another  line  is  printed.  When  the  group 
has  been  finished,  two  lines  are  printed  in  rapid  succession. 

Fission  Source  Program 

Stop  Printout 

52  "End  of  Run - "  Several  lines  of  information  are 

printed.  Unless  Edit  is  requested  by  depressing  switch  6 
or  the  iterations  have  converged  the  program  returns  to 


the  Group  Source  Program 
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Edit  Program 
Stop 
6253 

4323 

Other  Stops 
Stop 
161 

52 


2106 

2315 


Printout 

"Tape _ is  free.  Set  a  tape  10  for  output.  Depress 

start. " 

"End  —  Save  and  label  tape _ =  A  and  Tape _ 

=  E." 

Description 

Master  Tape  Generation . 

If  Vp  is  outside  tolerance,  follow  instructions  from  cust 
omer.  In  all  other  cases,  follow  instructions  on  the  printer. 
Group  Source  or  Group  Sweep. 

Follow  instructions  on  the  printer .  If  there  is  no  printout, 
record  index  register  C  and  push  start.  H  the  stop  occurs 
again,  display  index  register  C  and,  use  Emergency 
Restart . 

Group  Sweep. 

Use  Emergency  Restart. 

Group  Sweep. 

See  Restart  procedure. 
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Stop  Description 

4323  Edit, 

The  end.  Follow  instructions  on  the  printer, 

6253  Edit. 

Follow  instructions  on  the  printer . 

11646  Re-edit. 

Restart  Re-edit, 
any* others  Machine  error. 

If  a  stc^  occurs  before  the  Emergency  Restart  deck  has 
been  punched,  use  the  Control  Restart,  if  available;  other¬ 
wise,  start  TRIXY  over  from  the  beginning.  If  this  stop 
is  not  in  Edit,  use  the  Emergency  Restart.  If  this  stop 
is  in  Edit,  save  and  label  the  master  tape  and  tape  E  for 
Re-edit.  List  the  output  tape  (tape  10). 

UNUSUAL  OPERATING  PROCEDURE 


Loading  New  Liebmann  Factor 

1 .  Punch  card  for  new  oC  as  follows : 

.29,  ^  ,  (Control  punch  must  be  a  -  end  of  record,  floating 
point.  Do  not  use  an  end  of  file  card.) 

2.  Put  deck  for  MDL9C  followed  by  data  card  for  in  card  reader.  Do 
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not  put  BILD  on  front  of  MDL9C  deck.  The  MDL9C  deck  contains  its 
own  load  routine,  MBILD,  as  the  first  two  cards  of  the  deck.  Deck  is 
numbered  MDL9CN01  through  MDL9CN44. 

Caution;  There  are  three  different  pre- assembled  decks  for  MBILD  for 
the  "80C0",  ”16,000”  and  ”32,000”  word  memory.  Be  sure  to  use  the 
proper  deck.  The  identification  ”N”  punched  in  column  78  of  the  MDL9C 
deck  is  N=1  for  ”8000”,  N=2  for  ”16,000”,  and  N=3  for  ”32,000. 

3 .  Ready  card  reader . 

4.  Put  sense  switch  4  down. 

5.  After  data  have  been  loaded,  the  program  stops  at  (52)g.  Flow  instruct 
ions  on  printer.  Put  sense  switch  #4  up  before  pressing  start. 

RESTART  PROCEDURES 

A.  Restart  Procedure  (for  Group  Sweep). 

To  interrupt  TRIXY  in  the  middle  of  Group  Sweep: 

1.  Put  sense  switch  5  down. 

2.  When  stop  occurs  at  (2315)  ,  follow  instructions  on  printer  for 

O 

removing TRDCY  from  machine. 

NOTE  -  "Marking  each  tape  reel  by  number”  means  that  the  tape 
reel  on  unit  1  should  be  marked  1,  2,  3,  etc. 
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To  Restart  TRIXY; 

1 .  Replace  the  tape  reels  removed  from  the  machine  and  set  the 
numbers  on  the  units  to  those  indicated  by  the  markings  on  the  tape 
reel  (see  above). 

2.  Check  to  be  sure  step  1  has  been  done  properly. 

3.  If  this  is  the  first  time  TRIXY  has  been  restarted,  set  the  tape 
unit  that  is  now  tape  1  to  tape  3,  and  set  the  tape  unit  that  was 
tape  3  to  tape  1 . 

CAUTION  -  If  this  is  not  the  first  time  TRIXY  has  been  restarted, 
omit  this  step  and  leave  the  tape  numbers  as  they  are . 

4 .  Place  the  card  punched  out  when  TRIXY  was  interrupted,  the 
MDL9C  deck,  and  the  proper  Liebmann  factor  ( )  in  the  card 
reader .  The  value  for  must  be  read  from  the  printer  output 
of  the  Group  Sweep  Routine  (first  floating  point  number  printed 
on  a  line)  and  punched  on  a  special  card. 

5 .  Ready  the  card  reader . 

6.  Set  the  sense  switches  for  normal  operation  of  TRDCY  (see 

"Normal  Operating  Procedure”  p.  159  ). 

7.  Put  sense  switch  4  down  (  a  must  be  loaded). 


8. 


Press  "CLEAR"  and  "LOAD  TAPE", 
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B.  Emergency  Restart  (for  Group  Source  or  Sweep) 

This  restart  requires  that  the  master  tape  and  tape  E  be  saved. 
The  master  tape  is  on  unit  3  unless  the  Restart  described  in  Restart 
Procedure,  page  164,  has  been  used.  In  this  case  the  master  tape  is 
on  unit  1 .  Check  the  printer  to  see  if  a  run  has  been  completed.  If  a 
run  has  been  completed,  the  number  of  the  new  tape  E  is  indicated  in 
the  printout  from  the  end  of  the  last  run.  If  a  run  has  not  been  completed, 
tape  E  is  on  unit  5 . 

1 .  Set  the  master  tape  to  3  and  tape  E  to  5 . 

2.  Remove  the  "OMEGA”  deck  from  the  card  reader.  Take  one 
card  for  each  run  which  has  been  completed  from  the  front  of 
the  deck. 

3 .  Place  the  remainder  of  the  "OMEGA"  deck  in  back  of  the 
emergency  restart  deck. 

4 .  This  combined  deck  goes  into  the  card  reader . 

5.  Push  CLEAR  and  LOAD  CARDS. 

If  for  any  reason  it  is  impossible  to  continue,  save  and  label 
the  master  tape  and  tape  E ;  save  the  Emergency  Restart  and 
"OMEGA"  decks .  To  restart,  place  the  master  tape  and  tape  E 
on  the  machine  and  follow  instructions  1. - 5.  under 


Emergency  Restart. 
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C.  Control  Restart  (for  Master  Tape  Generation) 

This  restart  requires  that  tapes  on  units  1  and  2  be  saved. 

1.  Place  tapes  from  units  1  and  2  on  units  1  and  2  respectively. 

2.  Provide  blanks  on  units  3,  4,  5,  and  6. 

3.  Place  the  Control  Restart  deck  in  the  card  reader. 

4.  Push  CLEAR  and  LOAD  CARDS. 


